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Résumé 


Cette thèse s’intéresse à l’étude du nombre maximal de cycles limites des sys¬ 
tèmes différentiels ordinaires dépendant d’un petit paramètre. Plus particulièrement, 
on étudie deux classes de systèmes différentiels en utilisant la théorie de moyenni- 
sation. 

La première classe concerne l’étude des systèmes différentiels polynomiaux de 
Liénard généralisés de la forme : 

[x = y, 

\ÿ = -x - /(x) - g(x)y - h{x)y 2 - l(x)y 3 , 

où /(x) = efi(x) + £ 2 / 2 (x), g(x) = eg\{x) + £ 2 g 2 (x), h{x) = eh 1 (x) + £ 2 h 2 (x) et 
/(x) = eli(x) +£ 2 l 2 (x). Pour chaque k = 1,2, fk(x), gk(x), hk(x) et lk(x) sont de 
degrés n. e est un paramètre réel supposé petit. 

la deuxième classe concerne l’étude des systèmes différentiels de Kukles géné¬ 
ralisés de la forme : 

f x = —y + l(x), 

\ÿ = x - f(x) - g(x)y - h(x)y 2 - d 0 y 3 , 

où l(x) = sl 1 (x) + e 2 / 2 (x), /(x) = e/ 1 (x) + £ 2 f 2 (x), g(x) = £g 1 {x) + £ 2 g 2 (x), 
h(x) = f/i 1 (x) + £ 2 h 2 (x ) et do = £d\ + £ 2 c/q. Pour chaque k = 1,2, l k (x), / fc (x), 
g k (x) et h k \x) sont de degrés m, ni, n 2 et respectivement, ^ 0 est un nombre 
réel et £ est un petit paramètre. 

L’étude des deux classes est illustrée par des exemples. 

Mots clés : Cycle limite, théorie de moyennisation, système de Liénard, système 
de Kukles. 



Abstract 


This thesis is interested to the study of the maximum number of limit cycles of 
ordinary differential Systems depending of a small parameter. Using the averaging 
theory of first and second order, we study two classes of differential Systems. 

The first class deals with generalized polynomial Liénard differential System of 
the form 

\x = y , 

\y = -x - f(x) - g(x)y - h{x)y 2 - l(x)y 3 , 

where f(x) = efi(x) + £ 2 f 2 (x), g(x) = £gi(x) + £ 2 g 2 (x), h(x) = £hi(x) + £ 2 h 2 (x) 
and l(x) = £l\{x) + £ 2 1 2 {x) where fk{x), gk(x), hk(x ) and lk(x) hâve degree n for 
each k = 1,2, and e is a small parameter. 

The second class studied the generalized polynomial Kukles differential System of 
the form 

(x = -y + l(x ), 

\ÿ = X - f(x) - g(x)y - h(x)y 2 - d 0 y 3 , 

where l(x) = £Z 1 (x) + £ 2 l 2 (x), f(x) = e/ 1 (x) + £ 2 f 2 (x), g(x) = £g l {x) + £ 2 g 2 (x), 
h(x) = £h 1 (x) + £ 2 h 2 (x) and do = £d\ + £ 2 c(q where l k (x), f k (x), g k (x) and h k (x) 
hâve degree m, ni, n 2 and n 3 respectively, cIq 7^ 0 is a real number for each k = 1,2, 
and e is a small parameter. 

This study is illustrated by applications. 

Keywords : Limit cycle, Averaging theory, Liénard Systems, Kukles Systems. 
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Introduction générale 


Un système dynamique est un système physique qui évolue. Il peut évoluer 
dans le temps ou par rapport à une autre variable suivant l’espace de phases 
considéré. La trajectoire d’un objet en mouvement dans le temps est donc un 
système dynamique, ainsi que le nombre d’individu d’une population quelconque 
dans le temps, ou encore les valeurs d’une fonction (par exemple : y = 2x) 
par rapport à la valeur de x. Il est possible de suivre l’évolution de l’état d’un 
système physique dans le temps. Pour cela, on construit d’abord un modèle 
avec les lois physiques et les paramètres nécessaires et suffisantes pour caractéri¬ 
ser le système. Ce modèle est bien souvent constitué par des équations différentielles. 

Les équations différentielles sont apparues la première fois vers la fin du 17 eme 
siècle dans les travaux de Isaac Newton, Leibniz et Bernoulli. Elles se sont produites 
comme conséquence normale des efforts de ces grands savants d’appliquer les 
nouvelles idées du calcul à certains problèmes en mécanique. Plus tard la théorie 
d’intégration des équations différentielles a été développée par des analystes et des 
mécaniciens comme Lagrange, Poisson, Hamilton et Liouville aux 18 eme siècle et 
igeme gicles. Pendant plus de 300 ans, les équations différentielles ont servi l’outil 
essentiel pour d’écrire et analyser des problèmes dans beaucoup de disciplines 
scientifiques (Mécanique, Géométrie, Physique,•••). L’importance de cette notion a 
motivé des générations de mathématiciens et d’autres scientifiques pour développer 
des méthodes pour étudier les propriétés de leurs solutions. Par exemple, pour la 
mécanique non linéaire, on considère qu’elle fut fondée à la fin du 19 eme siècle 
par le mathématicien français Poincaré (Sur les courbes définies par des équa¬ 
tions différentielles, 1881-1886 ; Les méthodes nouvelles de la mécanique céleste, 
1892-1899). Il y a lieu de citer aussi le mathématicien russe Lyapunov, fondateur 
de la théorie de la stabilité (Le problème général de la stabilité du mouvement, 1892). 

Dans les années 1920, le statisticien Lotka et le mathématicien Volterra 
élaborent un modèle qui décrit la dynamique de systèmes écologiques où cohabitent 
un prédateur et sa proie. De façon emblématique, il est souvent appliqué aux 
lynx et aux lièvres des neiges dont les recensements précis ont été établis par la 
compagnie de la baie d’Hudson au 19 eme siècle. Ce modèle est constitué d’une 
équation différentielle qui traduit l’évolution des populations de chaque animal. 
Ainsi, quand les lynx sont nombreux, la population de lièvres décroît, ce qui 
entraîne la diminution du nombre de lynx, puis l’augmentation de celui de lièvres, 
libérés de la pression des prédateurs.... Ce type d’équations n’est pas réservé à la 
biologie et il est également pertinent notamment en physique : au début du 20 eme 
siècle, Van der Pol décrit de cette façon les oscillations d’un circuit électrique doté 
d’une lampe dont la résistance dépend de l’intensité du courant qui y passe. 
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Dans les deux exemples, à mesure que le temps croît, le comportement du sys¬ 
tème a tendance à devenir périodique. Graphiquement, les courbes qui représentent 
les solutions s’approchent d’un cycle : on parle de "cycle limite". 

L’acte de naissance des cycles limites est le Mémoire sur les courbes définies 
par une équation différentielle écrit par Poincaré en 1881 [42]. Un cycle limite 
est une orbite périodique isolée dans l’ensemble de toutes les orbites périodiques 
d’une équation différentielle. Rigoureusement, le modèle de Volterra n’admet 
un cycle limite que dans une version modifiée par l’écologiste canadien Holling, 
mais nous n’entrerons pas dans les détails. D’autres systèmes, plus complexes, 
peuvent avoir plusieurs cycles limites. Le comportement asymptotique sera encore 
périodique, mais la convergence aura heu vers tel ou tel cycle limite selon la po¬ 
sition initiale. Il est important de pouvoir déterminer le nombre de ces cycles limites. 

Cette question est au cœur de la deuxième partie du 16 eme problème de Hilbert 
[21] qui a été posé par David Hilbert au Paris conférence "congrès international des 
mathématiciens en 1900", en même temps que les autres 22 problèmes. On peut la 
formulée de la façon suivante : Quel est le nombre maximal des cycles limites qui 
peut avoir un système de degré n dans le plan? 

ix =P n (x,y), 

=Qn{x,y). 

Le problème des cycles limites n’est toujours pas résolu 112 ans après son énoncé! 
Dulac (1923) [13] proposa une démonstration assurant que H(n) est fini pour tout 
n. Mais cette démonstration comportait une erreur. La résolution du Dulac a été 
faite d’une façon indépendante par Ecalle, Martinet et Moussu en 1987, Ilyashenko 
en 1991 et par Ecalle en 1992. Cette résolution permet de montrer que H(n ) < oo. 
Petrovsky et Landis (1957) crurent trouver la valeur de H( 2) mais ils s’aperçurent 
d’une erreur dans leur propre démonstration Landis et Petrovski, 1967) avant 
que celui-ci ne soit infirmée par un contre exemple de Shi (1982) dans lequel un 
système quadratique a 4 cycles limites. Ainsi, si H(n) est un nombre fini pour tout 
n la seule chose que l’on sache est que H( 2) >4 et H( 3) > 11 (Jibin et Chunfu, 
1985 ; Zoladek, 1995). Christopher et lloyd (1995) ont donné une borne inférieur au 
nombre H(n) : H(n ) > 2 log(n). 

Les chercheurs considèrent deux classes spéciales des équations différentielles : 
des systèmes de Liénard généralisés et des systèmes de Kukles généralisés. 

Système de Liénard polynomial s’écrit : 

=V ’ (1) 
[y = -g(x) - f (x)y, 

où f(x) et g(x) sont deux polynômes de degré n et m respectivement. Ce système 
est très intéressant pour étudier les cycles limites des systèmes planaires, nous 
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avons une version simplifiée du problème de Hilbert dans [49] (1998). En 1977, Lins 
et al. [28] a déclaré la conjecture que si f(x) est de degré n > 1 alors le système 
(1) a au plus [^] cycles limites. Ils ont prouvé cette conjecture pour n = 1,2. La 
conjecture pour n = 3 a été prouvé par Li Chengzi et Llibre dans [25] (2012). Pour 
n > 5 la conjecture est pas vrai, voir De Maesschalck et Dumortier [12] (2011) et 
Dumortier et al. [16] (2007). Donc, il reste à savoir si la conjecture est vraie ou non 
pour n = 4. 

Il y a beaucoup de résultats concernant l’existence de cycles limites de faible 
amplitude pour (1), voir [36] (1988). On note H(m, n ) le nombre maximal de cycles 
limites qui peuvent exister simultanément pour le système (1). L’histoire commence 
avec Liénard (1928), un ingénieur français qui établit un théorème d’existence et 
d’unicité d’une solution périodique pour le système portant son nom (pour une 
revue complète des résultats antérieurs aux années soixante-dix, voir le livre de 
Sansone et Conti [48] (1964)). 

- En 1928, Liénard [27] montra que si ru. = 1 et si F(x) = f 0 ' f(s)ds est une 
fonction continue et impaire qui a une unique racine positive en x = a et qui 
est strictement croissante pour x > a, alors le système (1) possède un unique 
cycle limite. 

- En 1975, Rychkov [45] a prouvé que pour des polynômes F{x) impairs et de 
degré 5 et si m = 1, alors (1) n’a au plus que 2 cycles limites. 

- En 1977, Lins, De Melo et Pugh [28] ont prouvé que si m = 3 et n = 1, alors 
il n’y a au plus qu’un cycle limite. Ils ont de plus donné les conditions pour 
que ce cycle existe. Enfin ils ont conjecturé que si g(x) = x, il ne pouvait y 
avoir plus de E( m ^) cycles limites. 

- En 1983, Xianwu [53] a prouvé la conjecture pour certains cas où m = 4 et 
n = 1. 

Dumortier et al. dans [17] et [14] ont prouvé que 77(3,1) = 1. 

En 1997, Dumortier et Li [15] ont prouvé que 77(2,2) = 1. 

En 1998, Coppel [11] a prouvé que 77(2,1) = 1. 

Notons par 77(m, n ) le nombre maximum de cycles limites de faible amplitude du 
système (1). 

Blows, Lloyd [2] et Lynch ([37], [38]) ont donné les résultats suivants : 

- Si g est impair, alors 77(m, n) = [n/2]. 

- Si / est pair, alors H (m, n ) = n, quelque soit g. 

-Si / est impair, alors 77(m, 2n + 1) = [(m — 2) / 2] + n. 

- Si g(x) = x + g e (x), où g e est pair, alors T7(2m, 2) = m. 

En 1998, Gasull et Torregrosa [19] ont trouvé des bornes supérieurs pour 77(7,6), 
77(6,7) et 77(7,7). En 1999, Lynch et Cristopher [10] ont développé une nouvelle 
méthode algébrique qui permet de déterminer la constante de Liapunov du système 
(1). Ils ont donné les résultats suivants : 

H(m, 2) = [(2m + 1 )/3]. 

77(2, n) = [(2n + l)/3]. 
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Ê(m, 3) = 2[(3m + 2)/8] pour tout 1 < m < 50. 

H(3,n) = 2[(3n + 2)/8] pour tout 1 < n < 50. 

H(4:,k) = H(k, 4) pour k = 6, 7,8,9 et H( 5, 6) = H( 6, 5). 

En 2006, Yu et Han [54] ont obtenu des valeurs précises de H(m,n ) = H(n,m ) 
pour n = 4, m = 10,11,12,13, n = 5, m = 6, 7, 8, 9 et n = 6, m = 5, 6. 

En utilisant des techniques différentes (intégrales abéliennes, méthodes de la moyen- 
nisation, la méthode de Melnikov,...) on a les bornes inférieurs suivantes des nombres 
de Hilbert. On a H( 2) > 4 ([7], [50], [56]), H( 3) > 13 [24], H{ 4) > 22 [8], H{ 6) > 35 
[52], H(n) > kn 2 lnn [9] et H(n ) > 4(n + l) 2 (l.442695/n(n + 1) — g) + n — | [26]. 

En 2010, Llibre et al. [31] ont calculé le nombre maximal de cycles limites 
n) des systèmes de Liénard (1), en utilisant la théorie de moyennisation 
d’ordre k pour fc = 1,2,3. 

En 2012, Llibre et Valls [34] ont étudié par la méthode de moyennisation du premier 
et du deuxième ordre les systèmes suivants : 

x = y- gi{x) - fi(x)y, ÿ = -x - g 2 (x ) - f 2 {x)y, 

où gi(x), fi(x), g- 2 (x) et f 2 (x) sont de degrés k , k m et n, respectivement. 

En 2013, Llibre et Valls ont étudié par la méthode de moyennisation du troisième 
ordre les systèmes suivants (voir [32] et [33], respectivement) : 

X = y - gi(x), ÿ = -x - g 2 (x ) - f 2 (x)y, 
x = y- fi(x)y, ÿ = -x - g 2 {x) - f 2 (x)y , 

où gi(x), fi(x), g 2 (x) et f 2 {x) sont de degrés fc, k m and n, respectivement. 


Système de Kukles polynomial s’écrit : 

(x =-y + l(x), 

\ÿ =Q(x,y), 

où Q(x,y) et l(x) sont des polynômes avec des coefficients réels. Lorsque l(x) = 0 
et Q(x, y) est un polynôme de degré 3, le système (2) coïncident avec le système 
classique de Kukles : 


x = -y, 

ÿ = x + a 0 y + aiæ 2 + a 2 xy + a 3 y 2 + a 4 x 3 + ci 5 x 2 y + a ë xy 2 + ary 3 . 


( 3 ) 


Ce système a été introduit par Kukles dans [23] (1944) qui a donné des conditions 
nécessaires et suffisantes pour que le système (3) dispose d’un centre à l’origine. Ce 
système cubique sans le terme y 3 est appelé système de Kukles réduit. Sadovskii 
[46] (2003) a résolu le problème du centre-foyer de ce système avec la condition 
a 2 a 7 0 et a prouvé qu’il peut avoir 7 cycles limites. Rousseau et al. [44] (1997) 
apparaît une description des bifurcations locales des périodes critiques au voisinage 
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d’un centre non dégénéré des systèmes de Kukles réduits. Zang et al. [55] (2008) ont 
étudié le nombre et la distribution des cycles limites pour une classe des systèmes 
de Kukles réduits sous perturbation cubique. Chavarriga et al. [6] (2004) ont étudié 
le nombre maximal de cycles limites de faible amplitude qui peut coexister avec des 
courbes algébriques invariantes dans un système de Kukles. 

Plan de la thèse : 

Cette thèse se compose de quatre chapitres et une annexe. 

/ Chapitre 1 : le premier chapitre est un rappel de quelques notions générales que 
nous avons utilisé dans ce travail. Nous commençons par définir les systèmes 
dynamiques, les points d’équilibre et leur nature, la linéarisation au voisinage 
d’un point d’équilibre. Ensuite nous introduisons la notion d’un cycle limite 
et l’amplitude d’un cycle limite d’un système planaire. 

/ Chapitre 2 : dans ce chapitre, nous présentons des théorèmes importants de la 
méthode de moyennisation avec des exemples. On applique ces théorèmes pour 
la recherche des cycles limites des deux problèmes étudiés dans les chapitres 
trois et quatre. 

/ Chapitre 3 : ce chapitre concerne l’étude du nombre maximal des cycles li¬ 
mites qui bifurquent d’un centre linéaire perturbé par une classe généralisée 
de systèmes différentiels de Liénard de la forme : 

\x = y, 

\ÿ = -x- f(x) - g{x)y - h(x)y 2 - l(x)y 3 , 

où f[x), g(x), h(x) et l{x) sont des polynômes de degrés n. En utilisons la 
méthode de moyennisation du premier et deuxième ordre. 

Le résultat de ce chapitre est soumis pour publication dans le journal "Diffe- 
rential Equations and Dynamical Systems". 

/ Chapitre 4 : dans le dernier chapitre, nous donnons le nombre maximal de 
cycles limites de certains systèmes de Kukles généralisés de la forme : 

(x =-y + l(x), 

= X - f(x) - g(x)y - h(x)y 2 - d 0 y 3 , 

où l(x), /(x), g(x) et h(x) sont des polynômes de degrés m, ni, 712 et n 3 
respectivement, do 0 est un nombre réel. En utilisons la méthode de moyen¬ 
nisation du premier et deuxième ordre. 

Le résultat de ce chapitre a été publié dans le journal "Differential Equations 
and Dynamical Systems". 

Mellahi N, Boulfoul A and Makhlouf A. Maximum Number of Limit 
Cycles for Generalized KiLkles Polynomial Differential Systems. Differential 
Equations and Dynamical Systems, pp. 1-22, (2016). 
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initiales . 13 


Ce chapitre contient quelques notions générales et principales pour l’étude qua¬ 
litative des systèmes dynamiques et des équations différentielles ordinaires. 

1.1 Systèmes dynamiques 

Définition 1 Un système dynamique sur W 1 est une application U : R + x M n — > 
M n telle que 

1. U(.,x ) : M + — W 1 est continue. 

2. U(t ,.) : M + —> M n est continue. 

3. ?7(0, x ) = x. 

4- U(t + s,x ) = U(t, U(s,x)) pour t, s G M + et x G M n . 

Exemple 1.1.1 Soit le système 


x = Ax, x(0) = xo, (1-1) 

où A est une matrice constante, t G M + et x G M n . La solution de (1.1) est donnée 
par 


x(t) = e At xo. 















1.2. Flot d’une équation différentielle 
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Le système (1.1) engendre un système dynamiqiLe 

[/:K + xM n —X M n , 

U(t, x) = e At XQ. 

1.2 Flot d’une équation différentielle 

Définition 2 Soit le système non linéaire 


x = /(x) (1.2) 

avec la condition initiale x(0) = xq, xq £ E, E est un sous ensemble ouvert de M n et 
f € C 1 (E). Soit <p(t, xo) la solution de (1.2). L’ensemble des applications ipt définies 
par 


<Pt(xo) = <p(t,x 0 ) 

est appelé le flot de l’équation différentielle (1.2). 

Remarque 1.2.1 Le flot est dit autonome si f ne dépend pas explicitement du 
temps t, sinon il est dit non autonome. 

1.3 Points d’équilibre et linéarisation 

1.3.1 Points d’équilibre 

Définition 3 Le point xo G R n est appelé un point critique ou point d’équilibre du 
système (1.2) s’il vérifie 


f(x o) = 0. 


1.3.2 Linéarisation des systèmes 
Définition 4 Considérons le système (1.2). Le système 

x = Ax où A = Df(x 0 ) = ( fl 1 ' ] (1.3) 

\OXj J !<j 

est appelé Linéarisation de (1.2) en xq. 

Remarque 1.3.1 La linéarisation d’un système différentiel nous amène à l’étude 
de la nature des points critiques. 

Définition 5 Le point critique xq est dit hyperbolique si aucune des valeurs propres 
de la matrice jacobienne Df(x o) n’a de partie réelle nulle. 
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- 2 , xT(P~ - 1 , > 




Figure 1.1 - (0,0) est un point selle 



Figure 1.2 — (0,0) est un noeud stable 


1.3.3 Nature des points d’équilibre 

Soit le système différentiel linéaire (1.3) où A est une matrice 2 x 2 et soient Ai 
et A 2 les valeurs propres de cette matrice. On distingue les différents cas selon ces 
valeurs propres : 

1. Si Ai et A 2 sont réelles non milles et de signe différent, alors le point critique 
x = 0 est un point selle, il est toujours instable (voir FIGURE 1.1) 

2 . Si Ai et A 2 sont réelles de même signe, on a trois cas : 

(a) Si Ai < A 2 < 0, le point critique x = 0 est un nœud stable (voir FIGURE 

1 . 2 ). 

(b) Si 0 < Ai < A 2 , le point critique x = 0 est un nœud instable (voir FIGURE 

1.3). 

(c) Si Ai = A 2 = A 7 ^ 0, le point critique x = 0 est un nœud propre, il est 

stable si A < 0 et instable si A > 0 (voir FIGURE 1.4 et 1.5). 

3. Si Ai et A 2 sont des complexes conjuguées et Im {\\^) 7 ^ 0, alors le point 
critique x = 0 est un foyer. Il est stable si Re{X\p) < 0 et instable si Re( Ai^) > 
0 (voir FIGURE 1.6 et 1.7). 










1.3. Points d’équilibre et linéarisation 
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Figure 1.3 - (0,0) est un noeud instable 
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Figure 1.4 - (0,0) est un noeud propre stable 



Figure 1.5 - (0,0) est un noeud propre instable 
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Figure 1.7 - (0,0) est un foyer instable 











1.3. Points d’équilibre et linéarisation 
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Figure 1.8 - (0,0) est un centre 


4. Si Ai et À 2 sont imaginaires pures, alors le point critique x = 0 est un centre, 
il est stable mais pas asymptotiquement stable (voir FIGURE 1.8). 

1.3.4 Stabilité du point d’équilibre 

L’étude de la stabilité d’un point d’équilibre nous amène à connaître le compor¬ 
tement des trajectoires voisines de ce point d’équilibre. 

Définition 6 Soit le système 

dx 

— = f(t,x), x G M n , tel. (1.4) 

Supposons que f satisfait les conditions du théorème d’existence et d’unicité de la 
solution et soit 4>(t) la solution du système (1.4). On dit qu’un point d’équilibre p 
est stable si e > 0, 3 ô > 0 tel que si 

||$(to) ~P\\ < <5 =>- ||$(t) ~p\\ < à,Vt > t 0 . 

S’il existe de plus un voisinage de p tel que pour tout x dans ce voisinage lim <ï>(t) = 

t—¥ OO 

p alors le point d’équilibre p est dit asymptotiqiLement stable. 

Théorème 1.3.1 Soit xq un point d’équilibre pour le système (1.2) 

(i) Si toutes les valeurs propres de la matrice Jacobienne Df(x o) ont des parties 

réelles négatives, alors le point d’équilibre xq est asymptotiquement stable. 

( ii ) S’il existe au moins une valeur propre de Df(x o) avec une partie réelle positive, 

alors le point d’équilibre xq est instable. 

(iii) Si Df(xo) a des valeurs propres avec des parties réelles négatives et d’autres 
avec des parties réelles nulles, alors on ne peut rien dire siLr la stabilité du 
point d’équilibre xq. 
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1.4 Plan et portrait de phase 

Pour un système dont l’évolution au cours du temps t est décrite par la fonction 
à valeurs réelles x(t), on appelle trajectoire de phase une représentation géométrique 
cartésienne dans laquelle on reporte les positions au cours du temps t d’un point 
représentatif M d’abscisse x et d’ordonnée x = 

Cette terminologie est en accord avec celle de la physique statique : le plan ( x , x) 
ou le plan de phase s’identifie (avec p x = mx ) à l’espace de phase (x,p x ) du pro¬ 
blème, de telle sorte que la représentation au cours du temps du point M est bien 
la trajectoire du système dans son espace de phase. 

Une trajectoire de phase donnée est décrite à partir d’un point Mq(xq, xq) représen¬ 
tatif des conditions initiales de l’évolution considérée. L’ensemble des trajectoires 
de phase décrites par le système à partir de toutes les conditions initiales réalisables 
est le portrait de phase de celui-ci. 

Définition 7 Soit le système différentiel 

x = P{x,y), 

V = Q(x,y), 

où P et Q sont des polynômes en x et y à coefficients réels de degré d. 

Le portrait de phase est l’ensemble des orbites qui représentent les solutions du sys¬ 
tème (1.5) dans l’espace des phases ainsi que ces points critiques qui sont considérés 
comme des solutions constantes. Le plan (xoy) est appelé plan de phase. 


(1.5) 


1.5 Orbites périodiques et cycles limites 

1.5.1 Orbite périodique 

Définition 8 On appelle orbite périodique toute trajectoire <hj(x) de (1.5) telle qu’il 
existe un nombre T > 0, vérifiant 

!>(f + T, x) = <&(t,x). (1.6) 

Le plus petit réel T > 0 qui vérifie (1.6) est appelé période. 

1.5.2 Cycle limite 

Définition 9 Un cycle limite du système (1.5) est une orbite périodique isolée dans 
l’ensemble de toutes les orbites périodiques de ce système. 

Définition 10 L’amplitude du cycle limite est la valeur maximale de la variable x 
de ce cycle limite. 

Exemple 1.5.1 Soit le système 

J x = ax — y — ax(x 2 + y 2 ), 

\ y = x + ay- ay(x 2 + y 2 ), 


(1.7) 




1.6. Existence et unicité des solutions des problèmes à valeurs initialëtô 


tel que a est un paramètre. En coordonnées polaires x = r cos 6 et y = r sin0, le 
système (1.7) devient 

( r = ar{ 1 — r 2 ), 

\ 0 = 1 . 

Posons 


Alors 


f{r) = r = ar( 1 — r 2 ). 


f[r ) = 0 r = 0 ou r = 1. 

On a donc un point d’équilibre (0, 0) et un cycle limite d’amplitude r = 1 
(x(t),y(t)) = (cos (t + 0 O ), sin(t + 0 O )), 


où 


X 2 + y 2 = 1 . 


1. Pour a > 0, on a deux cas : 

(a) Si r < 1, alors f(r) > 0, d’où r{/*). 

(b) Si r > 1, alors f(r) < 0, d’où r(\). 

Donc le cycle limite d’amplitude r — 1 est stable. 

2. Pour a < 0 

(a) Si r < 1, alors f(r) < 0, d’où r(\). 

(b) St r > 1, alors f(r) > 0, d’où r{/ [ ). 

Donc le cycle limite d’amplitude r = 1 est instable. 

3. Si a = 0, le système a une infinité d’orbites périodiques, et il n’y a pas de cycle 
limite. 

Remarque 1.5.1 Les cycles limites apparaissent seulement dans les systèmes dif¬ 
férentiels non linéaires. 

1.6 Existence et unicité des solutions des problèmes à 
valeurs initiales 

Soit D un ouvert de M n , tç> < t < to + T, x G D et f(t , x) G M n . 
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Définition 11 On dit que f satisfait la condition de Lipschitz par rapport à x si 
||/(t,.Ti) - f(t,x 2 )Il < L\\xi - x 2 \\, 
où xi,x 2 £ D, et L est une constante. 

Nous pouvons maintenant formuler un théorème bien connu d’existence et d’unicité 
pour des problèmes à valeur initiale. 

Théorème 1.6.1 (existence et unicité) On considère le problème à valeur ini¬ 
tiale 


ilt = x(t 0 ) = x 0 , 

où x C D G R n , to < t < to + T. On suppose que : 

a) f(t, x) est continue par rapport à t , x sur G = [to, to + T) x D. 

b) f(t, x) satisfait la condition de Lipschitzienne en x. 

Alors le problème à valeur initiale admet une solution unique. 




Chapitre 2 


Méthode de moyennisation 


Sommaire 

2.1 Méthode de moyennisation du premier ordre. 16 

2.1.1 Application. 16 

2.2 Méthode de moyennisation du deuxième ordre. 19 

2.2.1 Application. 20 

2.3 Méthode de moyennisation du troisième ordre. 21 

2.3.1 Application. 22 

2.4 Autre méthode de moyennisation. 31 

2.4.1 Application. 32 


Dans ce chapitre, on s’intéresse à la méthode de moyennisation qui est appliquée 
aux systèmes de la forme : 

x = ef(x,t,e), (2.1) 

où t G I C R, x G M n , e un paramètre petit et / est T-périodique en t, l’équation 
moyennée associée à (2.1) s’écrit 

x = ( 2 . 2 ) 

où 

f°( x ) = -fJ o f(x,t,0)dt. (2.3) 

La recherche des racines positives du (2.3) réduit le problème de la détermination 

des solutions T-périodique de (2.1) qui est en général un problème difficile. 

La méthode de moyennisation est l’une des plus importantes méthodes per- 
turbatives utilisées actuellement dans l’étude des cycles limites des systèmes 
dynamiques. Cette méthode a une longue histoire qui commence avec les travaux 
classiques de Lagrange et Laplace en 1788 qui ont donné une justification intuitive 
de la méthode. La première formalisation de cette théorie a été faite par Fatou en 
1928. Elle a été introduite par Krylov et Bogoliubov en 1937 [22] et Bogoliubov 
et Mitropolskii [3] (1961). Elle a été ensuite développée par Verhulst [51] (1991), 
Marsden et McCracken [40] (1976), Sanders et Verhulst [47] (1985), Malkin [39] 
(1956) et Roseau [43] (1966), Llibre et Buica [4] (2004). D’autres formes et 
théorèmes de la méthode de moyennisation ont été démontrés ces dernières années 
et beaucoup d’articles ont été publiés concernant l’application de cette méthode 
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[5]. En utilisant la méthode de moyennisation, on réduit ce problème difficile des 
équations différentielles à la recherche des racines positives d’un système algébrique 
non linéaire. 

Dans ce chapitre, on présente des différents théorèmes de la théorie de moyen¬ 
nisation. 

2.1 Méthode de moyennisation du premier ordre 

Théorème 2.1.1 On considère le système suivant : 

_ dx _ x \ p x , £ \ (2.4) 

dt 

où Fi :KxD —> M n , R : M x Dx] — £/,£/[— > M n sont des fonctions continues, T- 
périodiques en la première variable et D un ouvert de M n . On définit F\o : D —> M n 
comme suit 

Fio(z) = j F 1 (s,z)ds, (2.5) 

supposons que 

(i) Fi et R sont localement lipschitziennes par rapport à x. 

(ii) Pour a G D tel que Fiq(o) = 0, il existe un voisinage V de a tel que Fifiz) fi- 0 
pour tout z G V\{a} et ds(Fio(a £ ), V, 0) fi 0. 

Alors, pour | e |> 0 suffisamment petit, il existe une solution du système (2.4) 

T-périodique isolée telle que —> a quand £ —> 0. 

Les hypothèses de ce théorème sont plus faibles que celui dans le théorème 11.5 de 
Verhulst [51], où à la place de (i) il suppose que : 

(j) Fi, R , D x Fi, D^Fi et D X R sont définies continues et bornées par une constante 

M (indépendante de £) dans [0, +oo [ X D, —£f < e < £f. 

A la place de (ii) il suppose que : 

(jj) pour a G D avec F\ifia) = 0, on a Jp 10 (a) fi 0, où D X F désigne la matrice 
jacobienne de F par rapport à x, D X F la matrice hessienne de F et «//(a) 
désigne le déterminant de la jacobienne de / calculée en a. 

2.1.1 Application 

1. L’équation différentielle de Van Der Pol : 

Soit l’équation de Van Der Pol 

x T x = £(1 — x 2 )x, (2.6) 


qui est équivalente au système différentiel 
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f x = y, 

\ ÿ = -X + e(l - x 2 )y. 


(2.7) 


En coordonnées polaires x; = rcos(0), y = r sin(0), le système perturbé (2.7), 
s’écrit sous la forme 


r = er sin 2 0(1 — r 2 cos 2 (0)), 

0 = — 1 + ecos0sin0(l — r 2 cos 2 0). 


D’ou 


On sait que 


D’ou 


dr er sin 2 0(1 — r 2 cos 2 0) 

d6 — 1 + e cos 0 sin 0(1 — r 2 cos 2 0) 


1 


1 — x 


1 + x + 0(x 2 ), 


x\ < 1. 


dr 

de 


= —er sin 2 0(1 — r 2 cos 2 0) + 0(e 2 ) 
= eF(d,r) + 0(e 2 ), 


De (2.5) on obtient 


1 f 2n 

F w (r) = — / F(6,r)d8 

^ J o 


( 2 . 8 ) 

(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 


Les cycles limite possibles pour l’équation (2.4) sont donnés par les racines 
positives de l’équation 

F w (r) = r -(r 2 -4) = 0. (2.12) 


Cette équation algébrique admet une seule solution positive r* = 2. Alors 
l’équation de Ven Der Pol (2.6) a un cycle limite d’amplitude r* = 2. 

1. Pour e > 0 ce cycle limite est stable car e^Fio(r*) < 0 (voir FIGURE 

2 . 1 ). 

2. Pour e < 0 ce cycle limite est instable car s-^Fioir*) > 0 (voir FIGURE 

2 . 2 ). 

2. Système différentielle de Liênard 

Soit le système de Liénard 


j x = y - e(aix + ... + a n x n ), 
\ ÿ = -x. 


(2.13) 
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Figure 2.1 - Cycle limite stable du système (2.7) pour e = 0.01. 



Figure 2.2 - Cycle limite instable du système (2.7) pour e = —0.01. 







2.2. Méthode de moyennisation du deuxième ordre 
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Le système (2.13) s’écrit en coordonnées polaires 

n 

r = —g akr k cos k+1 0, 
k= 1 

n 

6 = — 1 + g sin 6 ^ akr k ~ 1 cos k d. 


Ceci donne 


On prend 


k =1 


dr 

d9 


= £^2a k r k cos k+1 d + 0(e 2 ). 


k= 1 


-| /*27t i n p2n 

F io{r) = — F{6,r)d6 =-~y2a k r k cos k+1 8d6, 

27T J o 2îr Jo 


on aura 

, [^1 ,2. 

Fio(r) = — ^ a 2 fe+ir 2fe+1 / cos 2k+2 6dd. 

77 k =o 

Ce polynôme peut avoir au plus [-racines positives. Donc le système (2.13) 
peut avoir au plus cycles limites. 

2.2 Méthode de moyennisation du deuxième ordre 

Théorème 2.2.1 On considère le système différentiel 

x(t) = — = eFi(t, x) + g 2 F 2 (t, x ) + g 3 R(t, x, g), (2.14) 

dt 

où F\,F-2 :KxD —y R n , R : RxDx] — £f, £f [—> M n sont des fonctions continues, 
T-périodiques en la première variable et D un ouvert de M n . On définit F\q,F 2 q '■ 
D —> M n par (2.5) et 

F 2 o(z) = ^ J ^D x Fi(s,z) J Fi(t,z)dt + F 2 (s,z) S J ds, (2.15) 
supposons que 

(i) Pour tout t G M, F\ G C 1 , F \, _F 2 , R et D X F\ sont localement lipschitziennes 

par rapport à x. R est différentiable par rapport à g. 

(ii) pour V G D, un ensemble ouvert et borné et pour tout g G ] —g/, g/ [ 0, il existe 

a £ G V tel que F 10 (a £ ) + £F 20 {a £ ) = 0 et d B (F 10 (a £ ) + gF 20 (a £ , V, 0) f 0. 

Alors, pour |g| > 0 suffisamment petit, il existe une solution T-périodique isolée 
<£>(., g) de l’équation (2.14) telle que <p(0,£) = a £ . 
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Les conditions (i) et ( n ) de théorème 2.2.1 peuvent être remplacées par (j) et (jj) 
respectivement 

(j) Fi, R, D X F\, D 2 F\ et D X R sont définies continues et bornées par une constante 
M (indépendante de e) dans [0, +oo [ X D, —Ef < e < Ef. 


(jj) Fio(z) = 0, V 2; G D et pour a G D avec ^ 20 ( 0 ) = 0, on a Jp 20 (a) / 0. 


(2.16) 


2.2.1 Application 

Soit le système différentiel 

x = y + e(l — x 2 y ) + e 2 (æ — y 2 x), 
ÿ = —x + exy 2 + e 2 (l — y 2 x)y. 

En écrivant le système (2.16) en coordonnées polaires (r, #), on obtient 

f r = e(cosO + r 3 cos # sin 9 — 2r 3 cos 3 # sin #) + e 2 (r + 2r 4 cos 9 — r 3 cos 2 9 
—r 4 cos 5 9) — r 4 cos 0 + r 3 cos 4 9 ), 

0 = — 1 + e(—^1 + 2r 2 cos 2 0 — r 2 cos 4 0) + e 2 (—r 2 sin 9 cos 3 9 — r 3 cos 2 9 sin 9 


d’où 


ou 


+r 2 cos 9 sin 9 + r 3 sin 0 cos 0), 


d?’ 


Fi(9,r) 

F 2 (9,r) 


F- =EFi(9,r) + E 2 F 2 (9,r)FO(E :i ), 


2r 3 cos 3 d sin 9 — cos9 — r 3 cos 9 sin 9, 
r + 2r 4 cos 9 — r 3 cos 2 9 — r 4 cos 5 0) — r 4 cos 9 + r 3 cos 4 ( 

oin /9W o 1 y~i û I û „3 ^,-.,-,4 , 


+ COS - (2r 3 cos 2 d sin 9 — 1 — r 3 sin #)(— sin# + 2r 3 cos 2 9 — r 3 cos 4 d), 
r 

F 3 (9,r) = 0(e 3 ). 

On applique le théorème 2.2.1, donc nous avons 

1 f 2n 

Fio(r) = — / cos#(2?’ 3 cos 2 #sin# — 1 — r 3 sin#)d# 

27r do 


et 


= 0 , 

dFi 


dr 


(i 9,r) = 6r 2 cos 3 # sin# — 3r 2 cos# sin#, 


[0 1 1 

/ Fi(s,r)ds = —-?’ 3 cos 4 # — sin# + - cos 2 #). 

do 2 2 

Donc la fonction F 2 o(r) est 

-1 r 2 i r 0771 

d^oW = — y [- 7 ^(#, r) J F l (s,r)ds + F 2 (9,r)\d9 


= ~ 0 (r 2 ~ 8). 
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Figure 2.3 - Portrait de phase du système (2.16) pour e = 0.1. 


L’unique racine positive de F 2 q{t) = 0 est r = 2\/2. Le système différentiel (2.16) 
a un cycle limite instable d’amplitude r = 2\/2 car £F 2 0 (2\/2) = 2 > 0 
FIGURE 2.3 et FIGURE 2.4). 


voir 


2.3 Méthode de moyennisation du troisième ordre 


Théorème 2.3.1 On considère le système différentiel 

x(t) = = £ Fi( t,x) + e 2 F 2 {t, x) + e 3 F 3 (t,x) + e A FL{t, x, e), (2.17) 

où Fi,F 2 ,F 3 : K x P —y M. n ,R : M x Dx] — £/,£/[—> M n sont des fonctions 
continues, T-périodiques en la première variable et D un ouvert de M. n .On définit 
F 10 ,F 20 ,F 30 : D —> M n donnée par 


Fio(z) 

F2o(z) 

F 3 o(z) 


1 f 1 

T Jo Fl ( s ’ z ï ds ’ 


1 


D 


~Fi{s,z) i; F\(t, z)dt + F 2 (s , z)^j ds, 


\ D l F l{s,z) (yi(s, z)) 2 + D z Fi(s, z)y 2 (s, z) + F 3 (s, z) 


ds, 

(2.18) 


yi{s,z) = / Fi(t,z)dt, 

■J o 


y 2 (s,z) = / D z F 1 (t,z) / Fi(r,z)dr + F 2 (t, z) 

J o L J o 


dt, 


supposons que 
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Figure 2.4 - Portrait de phase du système (2.16) pour e = 0.01. 


(i) Fi(t ,.) G C 2 , F 2 (t ,.) eC 1 , Pour tout t £ R F \, F 2 , F 3 R , D X F\ et D X F 2 sont 

localement lipschitziennes par rapport à x. R est deux fois différentiable par 
rapport à e. 

(ii) pour V G D, un ensemble ouvert et borné et pour tout e G ] —£f,£f [ 0, il existe 
a £ eV tel que F 10 (a e ) + £F 20 (a £ ) + £ 2 F 30 (a £ ) = 0 et d B (F w (a £ ) + £F 20 (a £ + 
£ 2 Fs 0 (a £ ,V, 0) f 0. 

Alors, pour |e| > 0 suffsament petit, il existe une solution </?(.,£) T-périodique pour 
le système (2.17). 

Les conditions (i) et (ii) de théorème 2.3.2 peuvent être remplacées par (j) et (jj) 
respectivement 

(j) Fi, F 2 , R, D X F\, D 2 F\ et D X R sont définies continues et bornées par une 

constante M (indépendante de e) dans [0 ,+00 [ x D, —£f < £ < £f. 

(jj) Fio(z) = 0, F 2 q(z) = 0, V z G D et pour a G D avec F 3 q(ü) = 0, on a 
J F 30 ( a ) ^ 0- 


2.3.1 Application 

1. Cycles limites des équations différentiels polynomiaux de Liénard 
généralisés 

En 2009, Llibre, Mereu and Teixeira [31] ont étudié par la méthode de 
moyennisation, le nombre maximum de cycles limites H(m,n) qui peuvent 
bifurquer des orbites périodiques du centre linéaire x = y, ÿ = —x perturbé 
par une classe d’équations différentielles de Liénard de degré m et n de la 
forme : 
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x = y, 

ÿ = -X- Efc> 1 £ k {fn(x)y + Sm(*))> 


(2.19) 


où pour tout k le polynôme g^ n (x) est de degré m, le polynôme fh(x) est de 
degré n en x, £ est un petit paramètre. 

Ils ont calculé et donné des estimations inférieurs pour H(m,n). Plus pré¬ 
cisément, ils ont calculé le nombre maximum de cycles limites Hi i (m,n) qui 
peuvent bifurquer des orbites périodiques du centre linéaire x = y, ÿ = —x, 
en utilisant la méthode de moyennisation d’ordre k, pour k = 1,2,3. Ils ont 
obtenu aussi que Hk(m,n ) < H(m,n ) pour k = 1,2,3 pour les valeurs en 
lesquelles H ( m , n) est connu. 


Théorème 2.3.2 On suppose que pour k = 1,2,3 les polynômes fh(x) et 
sont de degré n et m respectivement où n,m > 1. Soit |e| > 0 est 
un paramètre suffisamment petit. Le nombre maximum de cycles limites des 
systèmes différentielles de Liénard généralisé (2.19) qui peuvent bifurquer du 
centre linéaire x = y, y = —x, en utilisant la théorème de moyennisation 

(a) du premier ordre est = [^]. 

(b) du deuxième ordre est É 2 (m,ri) = + [Ç], [§]}• 

(c) du troisième ordre est H^(m.,n) = [ n +™~ 1 ]. 

Preuve de la partie (a) du théorème 2.3.2 : Nous allons appliquer la mé¬ 
thode de moyennisation du premier ordre, pour cela nous écrivons le système 
(2.19) où k = 1, en coordonnées polaires (r, 6) où x = rcos#, y = rsin@. On 

n m 

pose /(x) = ^^OiX* et g(x) = ’^f iX 1 , Le système (2.19) s’écrit comme suit 
2=0 2=0 


r = —e 


Y. 

^i=0 


a, r î+1 cos* 6 sin 2 


6 + ^ bir 1 cos* 6 si 


55=0 


e = —i — 

r 


E 

\i =0 


a.;r* +1 cos î+1 


sm ( 


9 sin 9 + ^ bir * cos* +1 9 

i =0 / 


( 2 . 20 ) 


Si on prend 9 comme une nouvelle variable indépendante, le système (2.20) 
devient 


dr r 

d9 = 9 =£ 


a t cos* 9 sin 2 9 + ^ b {S cos* 9 sin 9 + 0(e 2 ), 


,!= 0 


i =0 


et 


r*27T / n 


Fio(r) = — 


2n 


yy a t cos* 9 sin 2 9 + bir 1 cos* 9 sin 9 d9. 


\i =0 


i =0 
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Pour calculer l’expression exacte de F\o(r). on utilise les formules suivantes 



r» 2n 


cos* 6 sin 9d9 


f a 2k , si i—2k, 

\ 0, si i=2k+l, 


0 , 


d’où 

[f] 

Fio(r) = £<&<*». ( 2 - 21 ) 

k=0 


Donc le polynôme Fio(r) a au plus [((] racines positives, et de plus on peut 
choisir les coefficients a 2 k de tel que F\o(r) admet exactement [^] racines 
positives simples. 


Preuve de la partie (b) du théorème 2.3.2 Nous allons appliquer 
la méthode de moyennisation du deuxième ordre. Pour cela nous écrivons 
le système (2.19) où k = 2, en coordonnées polaires (r, 9) où x = rcosd, 

n n n 

y = rsind, r > 0. On pose fi(x) = ^a;x*, f 2 (x) = ^qx*, g x (x) = ^TbjX 1 

2=0 2=0 2=0 

Tl 

et g(x) = ’y^djX 1 . 
i=0 

Le système (2.19) s’écrit comme suit 


n m 


E 


r= — s > a,r cos t/sm 


i =o 
f n 

2 / \ Qr^ 1 
1=0 


—e“ | ^2 Cir" ' " cos 


9 + ^2 bi rl cos* 9 sin 9 J 
i =o / 

m \ 

9 sin 2 9 + d.y cos* 9 sin 9 J , 


t=0 


9= 1 


1-( ^2 a i r%+1 cos* +1 9 sin 9 + ^ 6*r* cos* +1 9 


\ i =o 


i=0 


2 / n 

E 

i=0 


— — I ^2 CiX* +1 cos* +1 9 sin 9 + ^ dp’* cos* +1 9 


i=0 


( 2 . 22 ) 


Considérons 9 comme une nouvelle variable indépendante, alors le système 
(2.25) devient 
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ou 


Fi(r,6) 

fhM) 


E 

i =0 


air t+1 cos* 9 sin 2 9 + ^ ft^r* cos* 0 sin 9 , 


i=0 

m 


Cir* +1 cos* 0 sin 2 9 + ^ cifr* cos* 0 sin 0 


\i=0 


i=0 


—r cos 0 sin 9 ( ^ a;r* cos* 0 sin 9 + ^ 6jr* 1 cos* 9 J . 

\î=ü i =o / 

Déterminons la fonction F 2 o(r) correspondante. Pour cela, on pose Pio(r) = 0 
qui est L’équivalent au a 2 k = 0, ensuite calculons 


m m) 

dr 


y^(i + l)ajr* cos* 9 sin 2 


0 + ^ ï6jr* 1 


cos smt/, 


i=0 


i=0 


et 

yiO,#) 


= / F\(s,r)ds 

J o 

= air 2 (an sin# + «21 sin(30)) + ... + air l+l (au sin0 + a 2 ; sin(30) + ... 
+ a l± 3; sin((/ + 3)9) j + a o r(ai O 0 + «20 sin(2#)) + ... + a b r b+1 (ai b 9 

+a 2 &sin(2#) + ... + a&±4 6 sin((6 + 2)#)^ + b 0 ( 1 - cos0) + ... 
+b m r m (l — cos m+1 9), (2.23) 


où l est le plus grand nombre impair inférieur ou égal à n. b est le plus grand 
nombre pair inférieur ou égal à n. a 2 fc son t des constantes réelles qui figurent 
dans Jq 7 * cos 2fc 9 sin 2 9 pour tout k. 

On sait que d’après (2.21), Fio est identiquement nulle si et seulement si 
a 2 k = 0 pour tout k. On sait que : 


r* 27 T 


cos s 9 sin 2 9d9 


r* 27 T 


cos* 9 sin 9d9 


/ 0 

(*2tt 


cos* 9 sin 3 9d9 




cos 1 9 sirr 9 sin((2&; + 1 )9)d9 


r2i r 


cos* 5 #sin@sin((2/c + 1 )9)d9 


D’où 


( A 2i , si s=2i, 

\ 0, si s=2i+l. 

0. 

0. 

0. 

f fi 2 f +1 , si s=2i, 

\ 0, si s=2i+l. 



dF\(9,r) 

dr 


.y(9, r)d9 = r(âi 0 ai&o + (ài 2 ai& 2 + â 3 oa3&o)?’ 2 + 


••+ E 


cx.ij cii bj 
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n 

m 

l 

k 

i+/c-l 

2 

p+] + m 

impair 

pair 

n 

m 

n+m— 1 

2 

n—1 1 m 

2 ' 2 

pair 

pair 

n — 1 

m 

n— 1+ra— 1 

2 

((n-l)-l) . m 

2 T” 2 

impair 

impair 

n 

m — 1 

n+m—1—1 

2 

0-1) I {m- 1) 

2 1 2 

pair 

impair 

n — 1 

?7Z — 1 

n—1+m—1—1 

2 

((n-l)-l) (m-1) 

2 1 2 


Table 2.1 - Les valeures l + k — 1/2 ... 


où âij = -j±jA i+ j + i+j (auBj + ce 2 iBj + ... + pour tout i, j et 

k est le plus grand pair de n. 

De plus, on a 


r* 27 T 


F 2 (8,r)d8 = 


E 


j +i 


{•2tt 


dr ' / cos 

i=0 J ° 

i pair 

-2tt 


8 sin 2 8d8 + E E 2 r l+ i cubj 


j =0 2=1 
j pair i impair 


x / cos î+i+1 0 sin 2 8d8 = +)Cor + ... + +c&r 6+1 + 2(A 2 ai&cT 

J o 

++(<2360 + a\b 2 )r^ + ... + ++fc++ +fc E «i'À 

i+j=Z+fc 


On obtient alors le polynôme F 2 o(r) qui est égale à : 


r(p w aib 0 r + (pi 2 aii> 2 + /Osoû+o)?’ 2 + O 14 CH &4 + P 32 Û +2 + psoas&o)?’ 4 + ••• 
F pik<Hbkr l+k + AqCq + A 2 c 2 r~ + ... + +Q+), (2.24) 

où pij = àij + 2^4j_|_j_|_i pour tout i,j. Pour trouver les racines positives réelles 
de F 2 o(r), nous devons trouver les racines d’un polynôme dans r 2 de degré 
égal au max{++, §}. On a § = [§] et ++ = [^] + [f ] (voir TABLE 
2 . 1 ). 

Nous concluons que T 2 o a au plus max{[+l] + [y], [^]} racines positives. 

n 

Preuve de la partie (c) du théorème 2.3.2 On pose /i(x) = 

i= 0 

n n n n 

f2(x ) = = E Kæ *’ 5i(®) = E^+ = E d * æî et 

2=0 2=0 2=0 2=0 

n 

93 OO = 

i=0 

Le système (2.19) s’écrit comme suit 

J r = — sin# (eA + s 2 B + e 3 C) , 

[8 = -1- ( £ A + £ 2 B + £ 3 C) , 


(2.25) 
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où 


A = ^2 a?r* +1 cos* 0sin 2 0 + ^ bir 1 cos* #sin0, 

7=0 7=0 

n m 

B = ^ Cir l+l cos* 9 sin 2 9 + ^ di.r l cos* 9 sin 9, 

7=0 7=0 

n m 

C = ^ pir l+l cos* 0 sin 2 0 + ^ %r* cos* 9 sin 9. 

7=0 7=0 

Considérons 9 comme une nouvelle variable indépendante, alors le système 
(2.25) devient 


dr 

d9 


eA sin 9 + e 2 ( B sin 9 — 


A 2 cos 9 sin 9 


+ 


A 3 cos 2 9 sin 9 

2 2 AB cos 9 sin 9 


+ C sin 9 


De (2.25) Fio(r) est identiquement nulle si et seulement si ci 2 k = 0, pour tout 
k. De (2.24) nous obtenons i ? 2 o( r ) est identiquement nulle si et seulement si 
a.;, bj et Cfc qui satisfaisant l’équation suivante 


Cn — ^ 'y ^ Pijdi.bj, (2.26) 

^ i+j=fi +1 
i impair, j pair 


où p est pair. 

Pour appliquer la méthode de moyennissation du troisième ordre, nous devons 
calculer la fonction P 30 . Donc la preuve de la partie (c) du théorème 2.3.2 sera 
la conséquence directe des lennnes suivants. 

La preuve des lennnes suivants est simple et résulte de quelques calculs fasti¬ 
dieux. 


Lemme 2.3.1 Les fonctions correspondantes yi(9 , r) ety2(9,r) de la méthode 
de moyennisation du troisième ordre sont exprimées par (2.23) et 


2/2(0, r) = C 0 + C 1 r + C 2 r 2 + ... + C x r\ 




28 


Chapitre 2. Méthode de moyennisation 


Où À = max{2n + 1,2 m — 1} et 

C2k +1 = + X e° ai6 i é '+ X] fij a i a j@ 

i+j=2k i+j=2k+2 i+j=2k-\-l i+j=2k 

+^2fc+l + C2fc$ + S fe a 2m cos ( 2 * + 1 ) 6 ') + f 5Z 

i+jr=2fc+2 \i=0 / 


a,-a 




i î+jr=2fc 


/ fc +1 


y. WJ + y a i b j° +^ 2^+1 E a 2i+2 «+2? + 2 ) é ' 


i+.7=2fc+2 i+ < 7=2fc+l 

/ fc +1 


\ 2—0 


+ ^ dibj J y: a^ +1 sin(2z + 1)0 J + I ^ aibj + ^ aidjO 


i-\-j=2k-\-l \ i=0 


.i-\-j=2k+l i-\-j—2k 


+C 2 k) a\ i+ 2 sin(2i + 2)0^ , 


C 2 /,; = y c^OjOj + y djjbibj + y e^bjO [ y ciidj 


i-\-j=2k—l 


i-\-j=2k-\-\ 


i-\-j=2k 


yi-\-j=2k—l 


/ fc +1 


y 6++ y dibjo E b 2 i+i c°s(2 i + »» E bibj 


i+j=2k+l 
/ fc +1 


i-\-j=2k 


\i =0 


i i+j'— 2 fc+l 


E b 2i+2 cos ( 2i + 2)0 J + E dibj + C2k-1 + E dibj 9 


\i =0 
/fc+i 


^i+ji=2fc i+j=2k 

/ fc +1 


E ++1 sin(2i + 1)0 + y dibj y + +2 sin(2i + 2)0 


\i=0 / \ i J rj=2k 

J O Ü hl J J J fl 
i+ 

A 


\î=Ü 


Où a l 2i+l , a l 2i+2 , b 2i+1 ,b 2i+2 , c^-, <+, + , + sont des constantes réelles pour l = 
1,2 et k = 0,1,..., f. 

Lemme 2.3.2 L’intégrale Jq W r)(yi(s, r)) 2 ds est le polynôme 

7t(D 0 + D\r + D 2 r 2 + ... + D k r k ), (2.27) 


où 

si m>n+l et m ou n pair, 
si m>n+l et m ou n impair, 
si m <n + 1 et n pair, 
si m < n + 1 et n impair, 
et 

D x = y Pijk a i a j a k + y \jk a i b j b k + /* ] b ijk a i a j b ki 

i+j+k=x~ 1 i+j+fc=x+l î+j+fc=X 

pour y = 0,1,..., k où &lj k , \}j k , ôjj k sont des constantes réelles. 


k= ( 


n + 2m. — 1, 
n + 2m — 2, 
3n + 1, 

3 n, 
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Lemme 2.3.3 L’intégrale f^ 77 r)(y 2 (s, r)) 2 ds est le polynôme 


ou 


0 = < 


E 0 + E\r + E 2 r 2 + ... + 


( n + 2m, si m>n+l et n pair, 

n + 2 m — 1, si m>n+l et n impair, 

3n + 2, si m < n + 1 et n pair, 

( 3n + 1, si m < n + 1 et n impair, 


(2.28) 


et 


E-21+1 


'y ] Pijk a i a j a k + y ' \jk a i b j b k + 'y ] fiijbiCj 
i-\-j-\-k=2l —1 i-\-j -\-k = 2l-\-\ i-\-j=2l 

+ E + E u ïjk a i a J b kn, 

i-\-j=2l i+j+k=2l 


i pair 


E-21 — ^jf^ai.ajük + El \jk a i b j b k + El àij b i°j 

i-\-j-\-k=2l—2 i-\-j-\-k=2l i-\-j=2l—l 

+ E r iij a i d i+ E vïjk a i a jb k Tr + E QikUiCj'ïï, 

i+j=2l—l i-\-j-\-k=2l—l i-\-j+k=2l—2 

i pair i pair 

pour l = 0,1,..., | où \? jk , ôf jk , rjfj, v? jk , Q- k sont des constantes réelles. 
Lemme 2.3.4 L’intégrale f^ 77 ^^r{s,r)yi(s,r)ds est le polynôme 


|(F 0 + F\r + F 2 r 2 + ... + F v r v ), 

n + 2 m, si m>n+l et n pair, 

n + 2m — 1, si m>n+l et n impair, 

3 n + 2, si m < n + 1 et n pair, 

l, 3n + 1, si m < n + 1 et n impair, 


ou 


v = < 


et 

F 2 i +i 


(2.29) 


E Pijk a i a j a k + E ^jk a t b jh + E Ô L biC i 

i-\-j+k=2l—\ i-\-j+k=2l+l i+j=2l 

+ E vïjaidj, 


i+j=2l 


F-21 


E Pijk a i a j a k I E A |fc a * 6 A’ + E ô L biC i 


i-\-j-\-k=2l—2 i+j-\-k=2l i-\-j=2l—l 

+ E Vij a i d j "h E u ijk a i a j b k^ + E CfjfcOiCjTT, 

i+j=2Z—1 î+j+/c= 2Z— 1 i+j-|-fc=2Z—2 

i pair i pair 

pour l = 0,1,..., | où /3f jk , \f jk , uf jk , (f jk sont des constantes réelles. 
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Lemme 2.3.5 L’intégrale J Q 2?r Fo(s,r)ds est le polynôme 

— (Go + G\r + G 2 r 2 + ... + G^r^), (2.30) 

r 

où 

n + 2 m, si m>n+l et n pair, 

, n + 2m — 1 , si m>n+l et n impair, 

ip = 

3n + 2, si m < n + 1 et n pair, 

, 3n + l, si m <n + 1 et n impair, 

et 

°21 = $%k a i a i a k + X ijk a i b j b k + ÏMi 

i+j-\-k=2l—2 i-\-j+k=2l i-\-j=2l—l 

+ 22 VijOidj + p-21-2, 

i-\-j=2l—l 

pour l = 0,1,,.., | où / 3f jk , Xf jk , uf- k sont des constantes réelles. 

Par les lemmes 2.3.2, 2.3.3, 2.3.4 et 2.3.5 nous obtenons 


F 30 (r) = —( M 0 + M\r + M 2 r 2 + ... + M Q -\r Q 1 + M Q r e ), 


-1 l'2i~ i — ^ ' Btjk<Ha :! a k ^ ' Xijk^ibjbk 4“ ^ ' SijbîCj 

i+j-\-k=2l—l i-\-j+k=2l-\-l i-\-j=2l 

+ dijUidj + 22 v ij a i a j b k^, 

i-\-j=2l i+j=2l 

i pair 


M 2 I 


^ ^ fîijkQ'ibjbk ^ ^ ^ijk&i&j&k “1“ ^ ^ &ijCj 

i-\-j+k=2l i-\-j-\-k=2l—2 i-\-j=2l—l 

H - ^ ^ Tjij&idj + ^ ^ fJ j ijkQ'iQ j jQ j k “1“ ^2l—2P2l—2 

i+j=2l-l i+j+k=2l-2 

( \ 


+ 


^ ^ VijkQ j iQ j jbk H - ^ ^ Pijk&it 


i+j-\-k=2l—l 
K i pair 


i+j=2l-2 
i pair 


7T 


J 


+ £ 


T’ij k& i Q>j& k 'R i 


(2.31) 


i+ < 7+/c=2Z—2 
i pair 


pour l 

À 


= 0,1,..., | et 
n + 2m, 
n + 2m — 1, 
3n + 2, 

3n + 1, 


si m>n+l et n pair, 
si m>n+l et n impair, 
si m < n + 1 et n pair, 
si m < n + 1 et n impair. 
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On sait que a 2k = 0, on obtient Mq = 0 et M 2k+ \ = 0. De plus 
c fc — E i + j = k+i ciibj, pour tout k > b et M k = 0, pour tout k le plus grand 

i impair, j pair 

n + 2 m — 2, si nr>n+l et n, impair, 

n + 2m — 1, si m>n+l et n impair et n pair, 

3n + 2, si m < n + 1 et n,m pair, 

3n + 1, si m < n + 1 et n pair et m impair. 

Donc 


4o ir ) = ar(M 2 + Af 4 r 2 + M 6 r 4 + ... + Af A _ 4 r A 2 + Af A _ 2 r A ), 


où 

A4) = 22 P'ijk a ibj b k + 22 2jbiCj + 22 2ij a idj + 'CUlüPlj-2'K. 

i+j+k=u i+j=uj —1 i-\-j=uj—l 

i impair i pair i impair 

j pair j impair textjpair 

k impair 


Par conséquent, 40(4 es ^ un polynôme de degré À dans la variable r 2 . Donc 
F 30 (z) a au plus racines positives. 


2.4 Autre méthode de moyennisation 

On considère le problème de bifurcation des solutions T-périodiques du système 
différentiel 

x(t) = F 0 ( t , x) + eF\ ( t , x) + £ 2 F 2 ( t , x, e) , (2.32) 

avec e suffisamment petit. Les fonctions Fq, Fi : Ixfi -> R n , F 2 : MxOx (—eo, £o) ~^ 
R n sont de classe C 2 , T-périodique en t, et Q est un sous-ensenrble ouvert de M n . 
Le système non perturbé : 


x(t) = F 0 (t,x), (2.33) 

a une sous variété de dimension n des solutions périodiques. 

Soit x(t,z ) la solution du système non perturbé (2.33) telle que x(0,2:) = z. Nous 
écrivons la linéarisation du système non perturbé le long de la solution périodique 
x(t, z) est : 

ÿ = D X F 0 (t, x(t, z)) y. (2.34) 

Notons par M z ( t ) la matrice fondamentale du système différentiel linéaire (2.34) et 
par £ : W K X la projection de M n sur ses k premières coordonnées, c.à.d 

£ (xi, ...,x n ) = (xi, ...,x k ) ■ 

Une réponse au problème de la bifurcation des solutions périodiques pour le système 
(2.32) est donnée par les résultats suivants : 
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Théorème 2.4.1 Soit V C R k un ouvert borné et soit (3 : Cl(V ) —> M n k une 
fonction de classe C 2 . Supposons que : 

(i) Z = {z Q = (a, (3(a)),a G Cl(V)} C Q et pour chaque z a G Z la solution x(t, z a ) 
de (2.33) est T-périodique. 

{ii) Pour chaque z a G Z il y a une matrice fondamentale M z (t) de (2.34) telle 
que la matrice ALT 1 (0) — A/T 1 (27r) a dans le bloc supérieur droit matrice 
kx(n—k ) nulle et dans le bloc inférieur droit la matrice A a (n—k) x ( n—k)avec 
det(A Q ) / 0. 

On considère la fonction F : Cl(V) —> 

T'(a) = £ ^ J M" 1 (t) Fi {t, x{t, z a )) dt j . (2.35) 

S’il existe a G C avec F(a) = 0 et det {{dF/dz){a)) / 0, alors il existe une solution 
de période T <p{t, e) du système (2.32) tels que ip(t , e) —>• quand e —0. 

Nous supposons qu’il existe un ensemble ouvert V avec Cl{V) C fl tels que pour 
chaque z G Cl{V), x(t, z, 0) est T périodique. L’ensemble Cl(V ) est isochrone pour le 
système (2.32), c.à.d il est un ensemble formé seulement par des orbites périodiques, 
toutes ayant la même période. 

Corollaire 2.4.1 (perturbation d’un ensemble isochrone) un ensemble V ouvert et 
borné avec Cl(V ) C fl tel que pour chaque z G Cl(V), la solution x(t,z, 0) est T 
périodique, alors nous considérons la fonction T : Cl(V) —» M n 

T 

F(z) = j Mf 1 (t) Fi {t, x(t, z, 0)) dt. (2.36) 

o 

S’il existe a € V avec F(a) = 0 et det ((dF / dz)(a)) 0, alors il existe une solution 

T-périodique f(^f) du système (2.32) tels que (/?(0,e) —>■ a quand £ —> 0. 


2.4.1 Application 

1. Bifurcation de Hopf du système Michelson : 

En 2011, Llibre et Zhang [35] ont prouvé l’existence de bifurcation de zero-Hopf 
du système différentiel de Michelson 



(2.37) 


avec (x, y, z) G R 3 et c > 0 est un paramètre réel. Lorsque c = 0 le système 
(2.37) à l’origine comme une singularité ayant des valeurs propres 0, Fi. Et 
Lorsque c > 0 assez petit le système (2.37) a une orbite périodique qui tend 
vers l’origine quand c tend vers zéro. Maintenant nous présentons les résultats 


principaux. 
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Théorème 2.4.2 Poilt c > 0 suffisamment petit, le système de Michelson 
(2.37) a une bifurcation de zéro-hopf à l’origine pour c = 0. De plus pour 
c > 0 suffisamment petit, les cycles limites qui bifurquent satisfait x(t) = 
—2ccosf + O(c), y(t) = 2csinf + O(c) et z(t ) = 2ccotf + 0{c). 

Preuve : Pour £ ^ 0, on pose le changement des variables x = ex, y = eÿ, 
z = ez et c = ed, alors le système Michelson (2.37) devient : 

' x = y, 

< ÿ = z i (2.38) 

z = —y + ed 2 — e-x 2 . 
y 2 

Nous utilisons toujours x,y,z au lieu de x,ÿ,z. Posons x = x, y = rsin# et 
2 = r cos 9, alors le système (2.38) devient : 

x = r sin#, 

< r ="-{2d 2 - x 2 )cos9, (2.39) 

9 = 1- -(2d 2 -x 2 )sm6. 

2 r 

Ce système peut être écrit comme suit : 

rj np p 

— = r sin0 + -(2 d 2 — x 2 ) sin 2 9 + e 2 fi(9, r, e), 

{te 2 ( 2 ' 4 °) 

-æ 2 )cos 9 + E 2 f 2 {9,r,£). 

Où fi et f -2 sont des fonctions analytiques dans leurs variables. 

Pour (x'o. ro) / (0,0)) le système (2.40) £=0 a une solution 27T périodique. 

x(9) = ?’o + xq — ?’o cos 9 , r(9) = xq, (2-41) 


tel que x(0) = xq et r(0) = ro- La première équation variationnelle de (2.40) £=0 
le long de la solution (2.41) est 

/O sin 9 \ 

V 0 0 J \y 2 J ' 

Il a la matrice M(t) de solution fondamentale 



1 1 — cos 9 \ 

0 1 )’ 


(2.42) 


qui est indépendante de la condition initiale (xq, ro). En appliquant le corollaire 
2.4.1 au système différentiel (2.40), nous avons 


Hxo,ro) 


1 

2 



M~ l 


(2 d? - x 2 ) sin 2 9 
(2 d 2 — x 2 ) cos 9 


1(2.41 )d0. 
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Alors F(x 0 ,r 0 ) = (gi{x 0 ,r 0 ),g 2 (x 0 ,r 0 )) avec 

gi(x 0 ,r 0 ) = -{Ad 2 - 5 rl - 6r 0 z 0 - 2xq), gi{x 0 ,r 0 ) = ^r 0 (x 0 + r 0 ). 

Le système d’équations T = 0 a une seule solution réelle donnée par 
(z 0 ,r 0 ) = (—2d, 2d) et comme le jacobien DF{x Q , ro)|(x 0 ,r 0 )=(- 2 d, 2 d) = d 2 fi 0. 
Pour d > 0 et |e| >0 suffisamment petit et d’après le corollaire 2.4.1, le 
système (2.40) a une seule solution périodique (x(9,£),r(0,£)) de période 27T 
tel que (x(0,e), r(0, e)) — > {—2d,2d) quand e —> 0. les valeurs propres de 
£)J’(xo,ro)| (a . 0iro)=( _2d,2d) sont idi. Cela montre que l’orbite périodique est 
linéairement stable. 

Pour c > 0 suffisamment petit le système de Michelson (2.38) a une orbite 
périodique donnée par x{t) = —2c cos t + o(c), y(t) = 2c sin t + o(c) et 
z(t) = 2c cos t + o(c). 


2. Solution périodique pour une classe des équations différentielles non 
autonome de Newton : 

En 2016, Llibre et Makhlouf [29] ont étudié l’existence des solutions pério¬ 
diques de l’équation différentielle non autonome du deuxième ordre. 


x = —V x V{t, x), (2.43) 

Il II 2 

où V(t,x) = Ht— \-£W(t,x) avec W{t,x) une fonction 27T périodique par 
rapport à la variable t, £ est petit paramètre réel, x = (x\, ...,x n ) et 


V x V{x,t) 


( SV 

\ dx i : 


dV 

’’ dx„ 


Théorème 2.4.3 Nous définissons les fonctions 


d~k 


1 

2tt 


2t r 


/ 


dW 

dx k 


(' t , xio cos t + yio sin t, ...,x n0 cos t + y n0 sin t) dt, 


N-lk 


2tt 

i r dw 

~2nj COSt d^ k 


( t , Zio cos t + 2/10 sin t, ..., x n o cos t + y n o sinf) dt, 


o 

pour k = 1, ...,n. 5z la fonction W (t,x i, ...,x n ) est 2 it- périodique par rapport 
à la variable t, alors pour chaque £ fi t) suffisamment petit et pour chaqiLe 
(z* 0 ,..., x* 0 ,2/î 0 , •••, y^o) solution du système 


(æio, -, X n0 , 2/10, —,2/no) = 0, 

^ 2 fc (*^10, x n Q, 2 / 10 , •••, 2/no) = b, ^our k-l,...,n, 


(2.44) 
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satisfaisant 


det 


d(Ji ,..;fF2n) 


d (xio, -, x n o, yio, 


■,Uno) V^lO,-,a:nO,!/lO,-,yno)=(a:îo,-,<o,yî'o.-,<o)) ^ 0 ’ 


(2.45) 

l’équation différentielle non autonome (2.43) a une solution 2 tt périodique 
x (t, e) qui tend vers la solution 2 tt périodique (a^ 0 cos t + y\ 0 sin t, ..., x* 0 cos t 
+y* 0 sin t) de x + x = 0 quand e —>• 0. 


Preuve : Si y l = x*. pour i = 1,..., n, alors l’équation différentielle non auto¬ 
nome du deuxième ordre (2.43) peut être écrite comme un système différentiel 
de premier ordre suivant dans M 2n 


Xi = yu 

Vi = ~Xi - £ 


dW(t, x) 
dxi 


(2.46) 


Pour i = 1, ..., n, pour e = 0, (x, y) = (0.0) est le point d’équilibre de système 
(2.46). Les valeurs propres du système linéarisé on ce point sont toutes les 
imaginaires ±i,...,±i. La solution (x(t),y(t)) du système non perturbé (c.à.d 
système (2.46) avec s = 0) tels que (x(0),y(0)) = (xo,yo) es t 


Xi = x i0 cos t + y i0 sin t, 
Vi = Uio cos t — x i0 sin t. 


(2.47) 


Pour i = 1, ...n, toutes ces orbites sont 27T périodiques. On utilise la notation 
présentée dans la section 3, nous avons X = ( x,y),z = (xo, yo), Fo(x, t) = 
(y, — x), F\(x, t) = (0 ,—S7 x W(t,x)) et F 2 (x,t,e ) = (0,0). La solution de la 
matrice fondamentale (t) est indépendante de z et nous la dénoterons 
par M(t). Un calcul facile montre que : 


cos t 0 
0 cos t 


0 sint 0 ... 0 

0 0 sin t ... 0 


M(t) 


0 0 

— sin t 0 

0 — sin t 


cos t 0 0 ... sin t 

0 cos t 0 ... 0 

0 0 cos t ... 0 


— sin t 0 0 


0 


0 


cost 
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et 


cos t 0 
0 cos t 


0 — sin t 0 

0 0 — sin t 


0 

0 


M~\t) = 


0 0 
sin t 0 
0 sini 


cosf 

0 

0 


0 

cos t 
0 


0 

0 

cos t 


— sin t 
0 
0 


0 0 


sin t 0 0 


cos t 


Selon le corollaire 2.4.1 nous étudions les zéros a = (x'o, ijq) des composants 2 n 
de la fonction J-(a) donnée dans (3.17). Plus précisément nous avons J-(a) = 
(Fi(a), ...,F 2n (a)), tels que 


Fk{ot) 


F 2 k {a) 


1 f dW , 

— / sini-— {t,xiocost + yiosmt,...,x n ocost + y n osmt)dt, 
2tt J dx k 

o 

2tt 

1 f dW . . . ,, 

/ cos t— — (t, x w cos t + 2/10 suit, ...,x n0 cos t. + y n0 suit) dt, 
2vr J dx k 


pour k = 1,..., n. 

En illustrant par exemple : 

Considérons l’équation différentielle non autonome du Newton (2.43) 
avec le potentiel V (t,xi,x 2 ) = Xl \ + eW(t,xi,x 2 ), où W(t,x\,x 2 ) = 
(ax i + bx 2 + cxf + dx'ï,) sin t. Alors nous avons 


F\ (xiO) x 2 q , y 10 y y 20 ) 
F 2 (x lo ,x 2 o,yio,y2o) 
F3(x lo ,x 20 ,yio,y2o) 
F4(xio,x 2 o,yi 0 ,y2o) 


^(4a + 3c(x^ 0 + 3yf 0 )), 

g (4 b + 3d(x2 0 + 3//2 0 ));. 

3 

~^ cx 102 / 10 ) 

--dx 2 oy 2 o, 


Si ac < 0 et bd < 0, alors le système F\ = F 2 = F 3 = F 4 = 0 a 16 solutions 
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( æ ïo> æ 20 >2/îo> Vto) données par 


±2\ --,±2W--,0,0 


3c 


3 ci 



Nous calculons les déterminants correspondants aux solutions 
( æ 10> æ 20> Vwi Vw) 

(dÇFu ?\ , 

\d(xio, X20, 2/105 2/20) ( Xlo > :I:20 ’3 /10 ^ 20 ) = ( :l; îo> ;I: 20^l0^2o)y ’ 


pour ces quatre solutions (x^q, *20’ 2/îo’ 2 / 20 ) son t j^abcd, —-^abcd, —j^abcd et 
jçCibcd respectivement. Donc, l’équation différentielle (2.43) a quatre solutions 
27r-périodiques qui tendent vers les solutions périodiques de l’équation non 
perturbée. 




Chapitre 3 


Cycles limites d’une classe de 
systèmes différentiels généralisés 

de type Liénard 


Sommaire 

3.1 Présentation du problème et résultats principaux. 38 

3.2 Preuves des résultats. 39 

3.2.1 Preuve du théorème 3.1.1. 39 

3.2.2 Preuve du théorème 3.1.2. 40 


3.1 Présentation du problème et résultats principaux 

En 2014, Garcia et al. [18] ont étudié le nombre maximal de cycles limites de 
faible amplitude pour les systèmes différentiels plus généralisés de Liénard 

x = y, ÿ = -x - f(x) - g{x)y - h(x)y 2 , (3.1) 

où f(x), g(x), et h(x) sont de degrés n peut avoir [|] cycles limites en utilisant la 
théorie de moyennisation du premier ordre et n cycles limites en utilisant la théorie 
de moyennisation du deuxième ordre. 


Dans ce chapitre, en utilisant la théorie de moyennisation du premier et du 
deuxième ordre, nous étudions le nombre maximal de cycles limites qui bifurquent 
des orbites périodiques du centre linéaire d’une classe de systèmes différentiels de 
Liénard plus généralisée que (3.1) 


x = y , 

y = -x - f{x) - g(x)y - h(x)y 2 - l(x)y 3 , 


(3.2) 


où f(x) = efi{x) + £ 2 f 2 (x), g(x) = egi(x) + e 2 g 2 (x), h(x) = ehi(x) + £ 2 h 2 {x) et 
l(x) = eli(x) +e 2 l 2 (x). Pour chaque k = 1,2, fk(x), gk(x ), h^(x) et lk(x) sont de 
degrés n. £ est un paramètre réel supposé petit. 

Remarque : Si l{x) = h{x) = 0, le système (3.2) coïncide avec les systèmes diffé¬ 
rentiels généralisés de Liénard (1). 

Nos principaux résultats sont les deux théorèmes suivants : 
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Théorème 3.1.1 Pour |e| > 0 suffisamment petit, le nombre maximal de cycles 
limites du système différentiel (3.2) qui bifurquent du centre linéaire x = —y, y = x 
en utilisant la théorie de moyennisation du premier ordre est 

~ 77 

H 1 (n) = [-] + l. 

Théorème 3.1.2 Pour |e| > 0 suffisamment petit, le nombre maximal de cycles 
limites du système différentiel (3.2) qui bifurquent du centre linéaire x = —y, y = x 
en utilisant la théorie de moyennisation du deuxième ordre est 

H 2 {n ) = n + 1 . 

3.2 Preuves des résultats 

3.2.1 Preuve du théorème 3.1.1 

Pour appliquer la moyennisation d’ordre un, on écrit le système (3.2), en coor¬ 
données polaires 


x = rcos9, y = rsin9 avec k = 1. (3.3) 

Posons : 

n n n n 

fi( x ) = Y a i x \ 9i{x) = Y bi x \ hi(x) = Y °i xl et h{x) = Y dix\ (3.4) 

i=0 i=0 i= 0 i=0 

le système (3.2) devient 

n 

r= sYl i a i riR i( e ) + bir i+1 Ti{6) + c i r i+2 5 i (0) + d ir i+3 U t (9)) , 

i=0 

0 = -1 - - (ay cos i+1 6 + by +1 R i+ i(9) + ar l+2 T i+ i{9) 

i =o 

+ d i r i + 3 S i+1 (9)), 

où 

Rj(9) = cos J 0sin$, 

Tj (9) = cos J 0sin 2 9 = cosW — cos^ +2 9, 

Sj(Q) = cos- 7 9 sin 3 9 = cosW sin 9 — cos^ +2 9 sin 9, 

Ui{9) = cos- 7 9 sin 4 9 = cosW — 2cos^ +2 9 + cos^ +A 9. 

Considérons maintenant 9 comme nouvelle variable indépendante, le système (3.5) 
s’écrit sous la forme suivante : 


Fi(0,r) = Y, (a ynm + b,r‘ +1 T,(6) + c,r ,+2 S,(e ) + d,r‘ +3 U,(t))) . (3.6) 

i=0 


OÙ 
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Alors 


1 n 

*io0) = — ^2 (ayJiW + by + 1 ïi( 2n) + Cir l+2 Ji( 2tt) + d;P +3 À(27r)) , (3.7) 


i =0 


OU 


r2n p2n r2iv 

J k (2ir)= R k (9)d9,ï k (2ir) = T k (9) dO, J k (2ir) = / S k (d)d0, 


i o 

/'2'Tt 


10 


10 


4(2vr)= / *7*(0)d0. 

Jo 

Pour calculer l’expression exacte de Pio, nous utilisons les expressions des 
intégrales de l’annexe A. Nous obtenons 


F\o{r) = ^Y.[ h ÿ^y i + d 2M‘2yr 

' j7r i =0 


2i+2 


= r 


E 

i=0 


«i 


2 i + 1 (* + 1 )! 


i * r ' + d2i 2(ÏT2) r 


2i+2 


(3.8) 


ou 


«fc = 3.5. (2k - 1), a k+ 1 = (2/c + l)a fc . 


Donc le polynôme i ? io(r) a au plus [|] + 1 racines positives. 

Exemple 3.2.1 Dans cette section, nous montrons que la limite supérieure H] (n) 
est atteinte dans le cas où 1 < n < 5. Les calculs ont été vérifiés à l’aide de Maple. 
Nous notons Ffi 0 (r) la fonction F\q(t)/v correspondant, un calcul adéquat permet 
d’obtenir 


• F îo( r ) = l b 0 + î d 0 r 2 . 

• *10 (0 = *îo(0 = \ b o + (jjdo + f&2 y 2 + Jëd 2 r 4 . 

• *io (0 = *io( r ) = \ h o + (fdo + |6 2 )r 2 + (jgd 2 + j^hV 4 + jfg d 4 r 6 . 

Sz nous fixons Ùq = 2, b 2 = ^^,64 = 8 , do = — §, d 2 = 8 and d 4 = Ap. 
Aozzs avons 


*10 ( r ) = 1 — Qui a exactement un racine positive. 

* 10 ( r ) = *io( r ) = è(f ~ r 2 )(3 — r 2 ), gzzz a exactement deux racines positives. 
Ff 0 (r) = Ef 0 (r) = |(1 — r 2 )(2 — r 2 )(3 — r 2 ), qui a exactement trois racines positives. 

3.2.2 Preuve du théorème 3.1.2 

Pour démontrer le théorème 3.1.2, nous utilisons la méthode de moyennisation 
d’ordre deux. Tenant compte de (3.4) et en posant 

n n n n 

f- 2 (x) = Y, Pix\ g 2 (x) = Y Qix\ h 2 (x) = Y s %x l and l 2 ( x ) = J] Vix\ 

2=0 2=0 2=0 2=0 
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avec k = 1, 2. Le changement en coordonnées polaires ( r,9 ), transforme le système 
(3.3) en 


r = s Y {(Hr'Riid) + by+'T^) + c*r* +2 ^(0) + dy^U^O)) 

i=0 

n 

—e 2 Y ( Py R i(0 ) + + Si r i+2 Si(0) + Vi r i+3 U^O)) , 

z=0 

n 

o= -i- £ ~Y ( airi cosi+16 + by +1 R i+ i(o)+ c y +2 T i+1 (e)+dy +3 s i+ 1 ( 0 )) 

r i=0 

2 n 

Y O 5 * 7 "* cos?+1 0 + ®r î+1 lîj + i(é')r + Sir î+2 Tt+i(6») + i; i r* +3 5j + i(É')) . 

i=0 


Considérons maintenant 9 comme une variable indépendante, on obtient 

J = j = eF 1 (9,r)+e 2 F 2 (9,r) + 0(e 3 ), 
d9 9 


(3.9) 


(3.10) 


ou 

n 

Fi(0,r) = (apTR*(0) + \{9) + c i r i+2 S i {9)+dy +3 U i {9)), (3.11) 

i =0 
n 

F 2 (9, r) = Y ( PyRi(0 ) + qy +1 Ti(9) + S y +2 S t (9 ) + vy +3 Ui(9 )) 

i=0 

n 

- ~ Y ( atr'Riie) + by +l Ti{9) + C y +2 Si{9) + dy +3 U l (9)) 


i=0 


x ( a y cos i+1 9 + by +1 R l+1 (9) + C y +2 T i+1 {6) + dy +3 S i+1 (0 )) . 

(3.12) 


Déterminons la fonction F 2 q correspondante. Pour cela, posons F\q = 0 qui est 
équivalent à 

bo = d 2 [n ] = 0, 

< ‘ (3.13) 

b 2 i = 2^TÏ^2î-2, 1 F i < [|], 

Tout d’abord, en utilisant (3.13) et en remplaçant dans (3.11) on obtient 
n 1^1 

Fi(r,e) = ^( O yiî i (0) + cy +2 ^(0)) + Y {b 2i+1 r 2i+2 T 2i+1 (e) 

i=0 i=0 

[-] 

+ d 2i+l r 2i+A U 2i+1 {9)) + Y d 2 i- 2 T 2i+l (r 2i _2 - %^T 2i (0)\ . 

i =0 \ 1 / 
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Alors 


n— 1 1 
2 


dFl dr ^ = + (* + 2)cir t+1 Si(0)) + ^ ((2* + 2)b 2 i + ir 2t+1 

i=0 i =0 


2=0 

[f] 


x T 2i+1 (9) + (2* + 4)d 2i+1 r 2i+3 C/ 2m (0)) + ]T(2?: + l)d 2i _ 2 r 2i 


i=0 


x (r 2i _ 2 - j -T 2 i(o)). 


(3.14) 


En tenant compte du fait que 

1 / 2i - 2 \ „ 4i + 1 

2 2 * -2 V i - 1 


2i \ e 2i + 2 


2 2 *+ 2 (2i - 1) V * + 1 


2i + 2 


2 2i (2i - 1) V * 

ensuite, en utilisant les intégrales de l’annexe A, nous obtenons 

î+i 


0 = 0 , 


r «-i i 

2 


y(0,r) = ^ ^ajrVj(0) + Cjr î+2 Ji(0)^ + ^ ( 6 2i+ ir 2î+2 ^ y*,; sin(2Z + 1)0 


i=0 


i =0 


î+2 


1=0 
i +1 


+d 2i+ ir 2î+4 ^ 77; sin(2Z + 1)# J + X] d 2i-2r 2l+1 ^ À,; sin(2Z0), 

(3.15) 


1=0 


i =0 


i=0 


OU 


lu = 


li,l 


= < 


Ai = < 


7i,«-7i+i,J) 0 < l < i, 

^ 7&+14+1? ^ ^ A 1* 

7 i,i ~ 2 li+i,i + 7i+2,«, 0 < l < i, 

— 2 7i+i,i+i + 7i+ 2 ,i+i> ^ = * + 1) 

7î+2,i+2, l = i + 2. 

A-U - s±iA,i + f£?/W, 0 < l < i - 1, 

_4î+l fl. . I 2i+2 fl 7 — 2 

2 j_lrt,î t 2 i _ 1 p î +i iî , i — t, 

l = i + 1. 


2^+2 /O 


Maintenant, nous déterminons la fonction correspondante 

F 2 o(r) = F^(r) + Fi 0 (r), 


avec 


y 1 / \ _ 1 f 27r dFi(0,r) 

F2 » (r) “2 nj 0 

i r 2it 

F2o(r)= 2 nj 0 F ^°^ de - 

Dans les lemmes qui suivent, nous calculons les intégrales F^r) et E| 0 (r). 
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Lemme 3.2.1 L’intégrale p 2 o( r ) es t un polynôme en r donnée par : 

(k + l)ira k+s+ i 


[f] [^1 


k +1 


*20 W = S a 2s b 2k+l^ 


,2s+2k+l/ 


i r~" 1 1 L (2s^2%inC Sr 

s=0 k =0 \ 1=0 

k+2 


2 k+s (2s + l)(fc + s + 2)! ' 
3 (k + 2)na k+s+1 


+ a 2 ,d 2k+l r*-™+H2sJ2%,n C u - + 1)( * + s + 3 )!' 

(.k + l)(4fc + 10s + 15)7TQ!fc+s+l 


fc +1 


+ C 2 sb 2 k+\r 2s+ ~ ’ +3 ((2s + 2) ^2 ïk,l' ïïK s,l 

1=0 

k+2 

+ C 2 S d 2 k+ir 2s+2k+5 ((2s + 2) ^ ïk,l nK 2 s r 


1=0 


s+l 


2 k + s + 1 {2 s + l)(2s + 3 )(k + s + 3)b 

3 (k + 2)(4 k + 14s + 21)7 tq/ c _|_ s _|_i , 
” 2 fc+s+2 (2s + l)(2s + 3)(k + s + 4)!' 

3(2s + l)ita k+s 


+ +de À,,-r c tt + 2fc+ , +1 ; 2s : 1 ; ( ^: +2)! ) 

+ c 2 M i 2 - 2 +^m + 3 ) g à ,„ â- h + 2 t jig:^: +3 „ >) . 


(3.16) 


Preuve. D’après les expressions (3.14) et (3.15) nous obtenons 

n n / »2i r p 2n 

Fj 0 (r) = VV (majr^- 1 / R i (0)J j (9)d9 + ia iCj r i+j+1 R i (9)J j {9)d9 

i =o j=o ^ Jo Jo 

p2n p2n \ 

+ (i + 2)aa jr i+ i +1 S i (0)J j (0)d0 + (i + 2)c i c j r i+ i +3 S l {0)J J (0)d0 
Jo Jo y 

f n — 1 1 

ni 2 J / k +1 />27t 

+ ^2 ^2 [ a ib2k+ir l+ 2 k+ 1 {iy^^k,i / Ri( 9 )sm(( 2 l + l) 9 )d 9 + ( 2 k + 2 ) 

i =o fc=0 V i=o •'° 

p2n k-\-l „2 tt 

x / T 2 fc+ 1 (0)J i (0)d0) + Q 6 2 A;+ ir i+2fc+3 ((^ + 2) Vtm / $(0) sin((2Z + l)0)d0 

i=0 4o 

p2n k-\-2 ç2tx 

+ ( 2 fe + 2 ) / r 2 fe+ 1 ( 0 )^( 0 )d 0 ) + a i d 2 fe+ ir' i+ 2 fe+ 3 (ïV^M / ^( 0 ) sin( 2 Z + l) 0 d 0 

i=0 do 

p2n k-\-2 „2i r 

+ (2fc + 4) / Cl 2fc+ i(0) J<(0)d0) + Ci d 2fc+1 r i+2fc+5 ((i + 2) V / 5i(0) 

i=0 do 

rü=li [ ”-i i 

l' 2 'iv _ 1 2 I L 2 J 

x sin(2Z + l)9d9 + (2fc + 4) / U2k+i{9)Ji(9)d9) + E E« 2t+2 > 

J 0 „—n 7„ n 


s=0 /c=0 

S +1 /-27T 

X ( 6 2 s+ift 2 fc+ir 2s+2fc+3 V 7 s ,I / T 2 fc+i(/9) sin(2Z + l)0d0 + d 2 a +i& 2 fc+i»' 

z=o do 


2s+2fc+5 
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s+2 s+1 

E ~îs,l / T 2k+1 (0)sin(2l + 1 )0d9) + (2k + 4 )r 2s+2k+5 (b 2s+1 d 2k+1 £ % 

i=o 1=0 


C 2n 


U 2k+ i(0 ) sin(2Z + 1)0dO + d 2s+1 d 2k+1 r 


1=0 

s ~\~ e ^‘ /*27T 

2s+2k+7 y 2 ~ s j / u 2k+1 ( 0 ) 

1=0 Jo 

k-\-l «27T 


x sm 


' ° ^ ^ i -L /*Z7T 

m(2Z + l)0d0)) + EE aid 2k . 2 r i+2k (i^hi / W) sin(2/0)d0 + (2k + 1) 

s =0 k =0 ' '- 1 


;=i 

f 2tt or. I o /*2îr 

x (/ T 2fc _ 2 (0) - ——- / T 2fc (0)J i (0)d0)) + c i d 2fc _ 2 r î+2fc+2 


'o 


2fc- 1 


'o 


/c+1 /»27T 


2/c + 2 
2k- 

(*2tt 
/O 


/*Z7T /*Z7T 

x ((* + 2)V/3 W / S l (0)sm(2W)d0 + (2k + l)( T 2k - 2 (O)Ji(e)d0 

l=l Jo Jo 

r, <<27T \ Itl IV] / S + l 

T / T 2fc (0)ji(0)d0)) + V V r 2s+2fe+2 d 2s _ 2 6 2fe+1 ((2fc + 2)VÀ,, 

1 J s=0 k=0 \ 1=0 

r>2 7T /c+1 /*27 t 

x / T 2jfc+1 (0)sin(2Z0)d0 + (2s + l)(V 7 M / T 2s _ 2 (0) sin(2/ + l)0d0 

J 0 , n J 0 


2s + 2 
2s-lJ 0 

(*2tv 
10 


r2n 


s+1 


T 2s (0) sin((2Z + l)0)d0)) + d 2s _ 2 d 2fc+ ir 2s+2/c+4 ((2fc + 4)^+, 

Z=0 

^+2 /»27T 


f‘ Z7T 'J' 1 ^ /* Z7T 

x / C/ 2 fe+iW sin(2Z0)d0 + (2s + 1)(V+ M / T 2s _ 2 (0) sin(2Z + l)0d0 

Jo ;=0 Jo 

,2. x [fl [fl 

/ T 2s (0) sin(2Z + l)0)d0)) +VV r 2s+2fc+1 (2s + l)d 2s _ 2 d 2fc _ 2 

Jo / TTnl-k 


2s+ 2 
2s - 1 

fc+i 

52 ^ 

1=0 


C 2n 


2s+ 2 


s =0 fc=0 
2tt 


i( / T 2s _ 2 (0)sin(2/0)d0- 

0 2S — i J 0 


T 2s (0) sin(2Z0)d0^ 


J 


Pour simplifier l’expression du polynôme iT+r). nous utilisons les intégrales de 
l’appendice A. On a 


. / 0 27r R i (0)J j (0)d0 = / 0 27r R i (6)J j (0)d0 = f 2n Si(û)Jj(0)dû = 
f 2n Si(0)Jj(0)d0 = f 2n Ti(0 ) sin(2 1 + l)0d0 = J 2n U t (d) sin(2 1 + l)0d0 = 
f 277 Ti(0) sin(2 W)d0 = f 2n Ui(0) sin(2 W)d0 = 0. 

. Jq 77 Ti(6)Jj(0)dQ = 

2iqü+fc(27r), if i=2k, j=2s, 

2 X 2 (J2fc(2vr) -/ 2fc+2s+2 (27r)), if i=2k, j=2s+'l, 

_ 25Tî^2fc+2s+2(27r), if i=2k+l, j=2s, 

0, if i=2k+l, j=2s+l. 
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. f 0 2n T i (e)J J (6)d6 = 


(2s+l)(2s+3)^ 2fc (^ 7r )’ 

(2s+2)(2s+4) -^ 2 fe(27r) — ( 2 <s+2) ^2fe+2s+2(2vr) 

* + (2s+4) ^2fc+2s+4 (27F), 

~(2s+l) ^+23+2(271-) + (23+3)^+23+4(271-), 

. 0 ’ 


if i=2k, j=2s, 


if i=2k, j=2s+l, 
if i=2k+l, j^2s, 
if i=2k+l, j=2s+l. 


. C u 2k+1 (e)Me)dd 


2 s -|_2'^2/c+2s+2(2 / 7t) 5 i—2s, 

0, ifi=2s+l. 


. J 2n u 2k+1 (e)M6)de = 

{ ~ (2s+l ) -^2fc+2s+2 (2tt) + ( 2s +3) ^2fc+2s+4(27r), if i^2s, 

0, if i=2s+l. 


. sin(2i + l)êd« = | J 0 "’ * : l ‘ t |’ + J 

. Si(9) sin(2/ + 1 )0dt> = | ‘ : ^’ + J 

. t jt,(D Mm»={\ : 2 * + *• 

• IT Si(») sin(2 W)d0 = | ’ K,J: lzll + 1 ' 


Lemme 3.2.2 L’intégrale P| 0 (r) est un polynôme en r donnée par : 

[fl 

E 7Ta ^ „2s+l 

2-(»+l)! fe 


F») = 


ni r n— 1 i 

2 J L 2 


+ 


s=0 


2(s + 2) 


V2 s r 


,2s+3^ 


x r 


,2s+2fc+l lïOL s+k +1 


TTOs+fc+l ^ 2s+2fc+3 

■+*(« + * + 2)! + “ 2 * <Wir 


EE («2s^2fc+l 

s =0 fc =0 

37 TO' s _|_/-_^i 


2 s+fc+1 (s + + 3)! 

, „ u J 2 s +2 k+-i 37 ra s+fc+ i „ J ^ 2 s + 2 fc +5 157r 

+ c 2 sb 2 k+ir 2 s + k + 1 (s + k + 3)1 + c 2sd 2 k+ir -, 


a 


s-\-k+l 


2s+k+2 

_j 2s+2fc+l 3?r (fe + 1 ) a s+A: 

(s + fc + 4)! 2 fc+i 2 s ^2 2 s + fc ~ 1 (2s — l)(s + fe + 2)! 

_ j 2s+2fc+3 37T (3fc — 2g + 4)g s+ fc \ 

2fc+1 2s “ 2 2 s+fc (2s — l)(s + fe + 3)!/ ’ ^ 




46 


Chapitre 3. Cycles limites d’une classe de systèmes différentiels 

généralisés de type Liênard 


Preuve. En remplaçant l’expression (3.13) dans (3.12), nous avons 


F 2 (r, 9) = J2 {Pir l R,(9) + qi r i+1 T t (9) + s^S^O) + Vi r i+3 U 2i+1 (9 )) 

i =0 

1 / n 1^1 

- - 52( a ir iR i(0) + Cir i+2 Si(0)) + (b2k + ir 2k+ ' 2 T 2k+1 (9) 


, i =0 


k =0 


+ d 2k+1 r 2k+4 U 2k+1 (0)) + d 2k - 2 r 2i+l (T 2k _ 2 {6) 


k =0 


2k+ 2 
2k - 1 


T 2 *(0)) 


[T] 


x ^(ûir* cos î+1 (é») + Cir l+2 T i+1 {6)) + ^ (& 2 *H-i 7 ' 2fc+2 -R 2 fc+2(0) 


. i =o 


fc =0 


+ 2 d 2k - 2 r 2k+1 {R 2k -i - *^R 2k+1 (9)) 


k =0 


2k - 1 


Alors 


rL / rZn rZn rZn 

F 2 0 (r) = Yipy Ri(9)d9 + qi r i+1 T i (0)d0 + s i r i+2 Si{0)d3 + Vi r i+3 
i=o V «/o J o J o 


r>Zn \ / nZiz 

x / Uj(0)d0 ) — ( aja J T t+J '~ 1 / i2j(0) cos- 7 " 1-1 0d6 + a.jC.,r* +J+1 

J ° ) i =0 i=o V ■'o 

p2n p2n 

x / iîi( 0 )T i+ 1 ( 0 )d# + Qa i r i+J ' +1 / 5i(0)cos i+ 1 0d0 + Ci Ci r i+j+3 

JO JO 

r 21=11 

/*27T \ 21 L 2 J / f‘2lX 

x / 5,(0)r i+1 (0)d0 -VV o i 6 2 fc+ir i+2fc+1 ( / T 2fc+ i(0) cos i+1 

•'o / i=o fc=o V 

p2tt tâir p2n 

+ / R2k+2(0)R i (9)de) + a i d 2k+1 r i+2k+3 ( U 2k+1 (0)co S i+1 0d9 + S 2k+2 (9) 


'o 


/o 


'0 


f*27T 


f*27T 


x R i (9)d9) + c i b 2k+1 r ï+2lt+3 ( T 2k+1 (0)Ti(9)d0 + R 2k+2 (0)S l (9)d0) 

J 0 JO 

[211 

r2n r2n \ n ^ 2-1 

+ c i d 2 fe +ir i+2fe+5 ( / C/ 2 fe +1 (9)Ti(9)d9 + / -WW^(fl)^) TV (M 23-2 

•'o / i =0 3=0 

r2iT p2n 9 S 9 f 27r 

x r i+2s ( R i (9)R 2s - 1 (9)d9 + T 2s _ 2 (0) cos i+1 0d9 -Ri{9) 

J 0 J 0 Z8 — 1 7 0 

p2n p2n 

x R 2s+1 (6)d6+ T 2s (0)cos i+1 9d0)) + Ci d 2s - 2 r i+2s+2 ( S i (9)R 2s . l {0)d9 

J o Jo 
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+ 


x 

+ 

X 

X 

X 

X 


f^ T 2 S -2(mme- 

f n — 1 r n— 1 1 

2 L 2 


/*27T /*27T 

(/ s i (e)R 2s+ i(e)de + T 2s (d)T t (e)de) 

J 0 JO 

L_ 2~“ L_ 2~ J / /-2tT 

E E 6 2 S+ l^+l^ +2fc+3 / r 2s +l(0)iÎ2fc + 2(0)^ + d2 S +l62fc + ir 2S+2A:+5 

s=0 fc=o V ■'o 

p2n r2n 

/ U 2s+l (0)R 2k+2 (0)d9 + b 2s+l d 2k+1 r 2s+2k+5 T 2s+1 (0)S 2k+2 (0)d0 

J O JO 


[f] [V 1 


d 2 s+l,ld 2 fc+l r 

f»27T 

/O 


2s+2fc+7 


rzn \ . \ " J 

/ U 2s+ i(9)S 2k + 2 {9)d9 J — E E (^2A-2&2fc+i 

J ° J s= O fc=0 

r 2w 2 s -I- 2 

2s+2fc+2( / T 2s _ 2 (0)i? 2fc+2 (0)etë + / J^O^Tafc+i^cÆ - —E 

Jo «/ O Z 5 — 1 

p2ir p2n 

(/ r 2s (@)^ 2fc+2 (0)d0+ / i? 2s+1 (0)r 2fe+1 (0)d0)) + c/ 2s _ 2 d 2fe+1 r 2s+2fc+4 

JO JO 

/•2-7I- /•2tt 2s + 2 r 27r 

(/ T 2s _ 2 (0)S 2k+2 (0)d0+ R 2s _ 1 (e)U 2k+1 (0)d0 - ~—r( T 2s (9 ) 

JO Jo JS — I Jo 


/•27T \ If [fl 

5 , 2fc+2(Ê , )dé ) + / i? 2s+ l W^2 fe+ 1 E E E ^-2^_ 2 r 

■'O / i=0 fc=0 


2 s — 1 ,/q 

2a+2/c+3 


/ r 27r 2 s 4- 2 / ,2it 

x n T 2s _ 2 (9)R 2k _ 1 (9)d9 - ^E(y o T 2s _ 2 (0).R 2fc+1 (0)dl9 

+ T 2s (0)R 2k _ 1 (9)d0) + |^^ 2 Çj T 2s (0)R 2k+l (9)d9^j . 

Pour une expression explicite du polynôme Pi^r), en utilisant les intégrales de 
l’appendice A nous avons 

• / 0 2w Æi(0)d# = / 27r Si(éOdé> = J 0 2 ^ cos^ 9Rj(9)d9 = f 2 * cos i 9Sj{0)d0 = 
fo” R i (9)T j (9)d9 = f 2n S i (0)T J (0)d0 = J 27T R i (9)U j (0)d9 = 
f 0 2w S i (0)U j (0)d0 = 0. 


. f o 2n T t (0)d9 


0, if i = 2/o + l, 

12,(2.) = 4^, ifi = 2fc. 


. J 2n cos i+1 0T 2fc+1 (tf)d0 = / 0 2?r = 

J 0 , if * = 2s + 1 , 

| J 2fc+2s+2 (27r) = 2t+ r^ +2) , , if i = 2s. 
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• C R2k + 2(0)Si(0)d0 = lf T 2k+l {8)T l+1 {8)d8 = f 2n cos i+1 9U 2k+1 (8)d8 = 
C S 2k+2 (8)Ri(8)dd 

( 0, if i = 2s + 1, 

{ ï 2k+2s+2 {2Tï) = 2fc+ 3 s +2(fc+t+3)! » if * = 2s - 


C U 2k+1 (8)T.i(8)d8 = f 2n S 2k+l (8)Si{8)d8 


), if i = 2s, 

2 fc+ 2 s+ 2 ( 2 vr) — ï 2k+2s+4 (2ir) = 2 fc+s+3(^ ) _^ + 4); i if i = 2s + 1. 


D’après les lemmes 3.2.1 et 3.2.2, nous avons 
F2o(r) = F^{r) + FUr) 


L j 

£ 


fc+i 


g 2s r 2s+1 + 


^2s+3 


rni r n — 1 

L 2 J L 2 


+ ££ («2s&2fc+l 


s =0 fc =0 


2(^T 2sr 

^ I 1 

v „2s+2k+l/ c ST' ~ f i _ (2s + k + 2)a k+s+ i | ^ 2s+2fc+3 

( 2^7 k,l S,l 2 fc + s+ 1 ( 2 s +l)(fc + s + 2 )!^ 2s 2k+1 


_ (12s 2 + Ak 2 + 34s + lOsk + 19k + 24)afc +s+1 

X S ’ 1 2 fc+s+2 (2s + l)(2s + 3 )(k + s + 3)! ’ 

|_2 

+ £ fcA, - > + 

s , 2s+2k+bt/ , -n = T s 3(20s 2 + 4fe 2 + 68s + 14sA; + 29fc + 57)afe +s+ i 

X " ((S + 1 ] ^ k ’ 1 °’ 1 -2 fc+s + 3 ( 2s + l)(2s + 3)(fc + s + 4)!- } 


r 2 s+ 2 k+i(( 2s + 1 ) 3 /=< | 3(4A: + 2s + 5)a fc + s x 

+ a 2k+1 d 2s - 2 i { ——^/3 s j Cktl + 2 fe+s+ 2( 2s _ i) (fc + s + 2)! } 


i=i 

, J 2s+2fc+3 / (^2 + 3 ) r if 3{12k + 2s + 2l)a k+s \ 

+ 2fc+i 2s—2 ( 2 2LjP s ,i k,l 2 k + s+3 (2s-l)(k + s + 3)VJ ' 


i=i 


Pour trouver les racines positives réelles de F 2 o(r), nous devons trouver les racines 
d’un polynôme dans r 2 de degré égal au [|] + f 21 ^] + 2. Puisque + [§] + 2 

prend la valeur f 22 ^] + f 22 ^] + 2 si n est impair et la valeur [|] + f 22 ^] + 2 si n 
est pair, soit n + 1 dans les deux cas. Nous concluons que F 2 o(r) a au plus n + 1 
racines positives simples. 
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Exemple 3.2.2 Dans cette section, nous montrons que la limite supérieure H 2 (n) 
est atteinte dans le cas où 1 < n < 5. Les calculs ont été vérifiés à l’aide de Maple. 
Nous notons F^r) la fonction F 2 o(r)/r correspondant, un calcul adéquat permet 
d’obtenir 

• F 2 o( r ) = 2 ® - la 0 h + (§v 0 - |c 0 6i - |dia 0 )r 2 - ^dic 0 r 4 . 

• F 2 o( r ) = ho ~ hoh + (î v o ~ - |dia 0 - \b\a 2 + \d 0 a 3 + \q 2 )r 2 + 

(-^c?ic 0 - ^bic 2 + J^doci + F v 2 - jQdia 2 )r A - î |g(iic 2 r 6 . 

• F $o( r ) = ho - 2°0^1 + (1^0 - |c 0 6 i - fdiao - |bia 2 + |d 0 ai + \q 2 - |a 0 b 3 )r 2 + 
(-JèdiGo - c 2 + Ygdoci + jqV 2 - pdia 2 - ^ c 0 b 3 + j^d 0 a 3 - j^a 0 d 3 - F b 3 a 2 )r 4 
+ (-,2 8 dic 2 - 4sb 3 C2 - m d 3 a 2 - m d 3 C 0 + m d 0 c 3)r 6 - îio d 3C2r 8 . 

• F A 0 (r ) = \q 0 - hobi + (f^o - |c 0 61 - §c?ia 0 - |&ia 2 + |d 0 ai + |g 2 - 
la 0 b 3 )r 2 + (^<74 - ^a 4 &i - ^dic 0 - ^&ic 2 + F d 0 c 1 + ^w 2 - F dia 2 - ^c 0 & 3 + 

45^0^3 Yé^ , od 3 jg&3®2 “1“ ygd 2 ai)r “I - ( 42 g94 H“ l 2 ga 3 d 2 128^3^1 128^ 4 ^ 4 

^250463 + Y 2 g(i 2 Ci — 72g d 1 C‘2 — j 2 gb 3 C 2 — Ï 2 sd 3 ü 2 — 22 gd 3 Co + 728doC3)r 6 + (— Ï 2 Ô()d 3 C 2 ~ 
42gQ C 4 (/i — 7280^4^3 + 256^2C 3 — Ï2gô a 4rf 3 )^ 8 — 1(524 C 4 d 3 r 10 • 

• F 20 (t ) = g9o - 2 a o^i + ( 1^0 - h°bi - i^i°o - I& 1 Û 2 + fdoai + §92 - 
hob 3 ) r2 + ( 7^4 - 75(1461 - ^dic 0 - ^&ic 2 + ^doci + ^2 - ^dia 2 - ^ c 0 b 3 + 
Tëd 0 a 3 - jQÜ 0 d 3 - F b 3 a 2 + j^d 2 ai - ^a 0 6 5 )r 4 + (7^4 + 7§ga 3 d2 - 7^461 - 

128 72 g 04^3 “ 1 “ 728^201 128 d4 ^2 128 ^3 ^'2 128 d 0^2 128 d 0^0 “ 1 “ 128 d 0^3 “I - 

Ï28^0 a 5 — 3g4&5°2 — 728 a 4^5 — ggjCo^)?’ + ( — 728Ü^3 C 2 — 7280 C4 ^ 1 — Î28Ü C4 ^ 3 + 

25gd 2 C 3 — 7280^4^3 — 7 ô 24&5 c 2 — 256 C 0^5 — 256^5 a 2 + 256 a 5 ^2 — 728Ü a 4^5) r8 + 

(3Ü72 d 2 C 5 ~ 70M C4 ^ 5 “ 75§ggC2 <^5 — 57§o a4 ^5 ~ TQ^^s)?’ 10 ~ 77M0 C4 ^5^ 12 - 

Si nous fixons qo = l,ç 2 = —|,&i = — 1 ,63 = |,ao = l,ai = 0,a 2 = — 2 , a 3 = 
-|îi,co = 2, Cl = l,c 2 = -77, C 3 = ^L,d 0 = — di = -|,d2 = l,d 3 = = -6 

et V 2 = Tp. ./Vous avons 

F 2 q( r) = |(1 — ? ,2 )(2 — r 2 ), qui a exactement deux racines positives. 

^ 2 o( r ) = g(l — r 2 )(2 — r 2 )(3 — r 2 ), qui a exactement trois racines positives. 

F 2 o( r ) = ^(l -r2 )(2 — r 2 )(3 — r 2 )(4 — r 2 ), qui a exactement quatre racines positives. 

D’autre part, si nous choisissons qo 
| 5 ao = l,ai = 0,a 2 = —2,a 3 = 

16 _ 244 „ _ 1216 j _ 8 j 

3 , C 3 — ^ 57 , C 4 — 375 . «0 — 3 , ai 

U 4 . = 40398 ^ 5 ^ • -/Voas avons 

F 2 o( r ) = î 5 ô (1 — r 2 )(2 — r 2 )(3 — r 2 )(4 — r 2 )(5 — r 2 ), qui a exactement cinq 
racines positives. 


1,92 = 

119 _ 

342 ’ “ 4 " 

= -s.* = 


34 „ _ 1537022 t _ _i h _ 
5,94 — H8125 — 1,03 — 

749522 „ _ o „ _ t „ _ 
118125’ C ° — Cl — 4 ’ C2 — 
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généralisés de type Liênard 


Enfin si nous choisissons go = l,g 2 = — ^,g 4 = yy, b\ = — 1,63 

1 ^_ o _ 419 _ ic _ _ 96528128 _ o _ 

l,ftl U, Û 2 A ^3 342 ? ^4 40,^5 520695 ’ ^0 A ^1 


_ l^lU _ I ^ ,y «J -J 

TÿT ’ L 1 — -g- , C5 — 2565 ’ uo ~ — 3 ’ U1 — 
V2 = yp et V4 = 15 g94g34 39 • Nous avons 

F 2 o ( r ) = 750 0 “ r 2 K 2 - r2 )( 3 “ r2 )( 4 - 
racines positives. 


— c?2 = 1,^3 = X5) cfe 

r 2 )(5 — ? ,2 )(6 — r 2 ), qui 


— §> &5 — 2 , ao — 
1,C 2 = -^,C 3 = 
= — = — ® 

a exactement six 
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4.1 Introduction et résultat principaux 


En 2011, Llibre et Mereu [30] ont étudié par la méthode de la moyennisation le 
nombre des cycles limites qui bifurquent des orbites périodiques d’un centre linéaire 
x = y, ÿ = —x, du système différentiel polynomial de Kukles de type : 


x = y, (4.1) 

y = -x - /(x) - g[x)y - h{x)y 2 - d 0 y 3 , 


où les polynômes /(x), g(x) et h{x) sont de degrés ni, ri 2 et n 3 respectivement, 
do 7 ^ 0 est un nombre réel. Ils ont prouvé que le système (4.1) peut avoir = max 
l} cycles limites en utilisant la théorie de moyennisation du premier ordre et 

(IH + [=¥3. Pf ]. m. PA + [¥]•[?]. PA] + [fl • 

[n 3 _KL] _|_ [nij ^ | j C y C l e s limites en utilisant la théorie de moyennisation du 
deuxième ordre. 


Dans ce chapitre, nous nous intéressons au nombre maximum des cycles limites 
d’une classe de systèmes différentiels de Kukles plus généralisée que (4.1) : 

x = 


y 


—y + l(x), 

x ~ f(x) - g(x)y - h(x)y 2 - d 0 y 3 


(4.2) 
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où l(x) = el 1 (x) + e 2 l 2 (x), f{x) = e/ 1 (x) + e 2 f 2 (x), g(x) = eg 1 {x) + s 2 g 2 (x), 
h(x) = eh 1 (x) + e 2 h?(x) et do = ed g + e 2 cig. Pour chaque k = 1,2, l k (x), f k (x), 
g k (x) et h k \x) sont de degrés m, m, n 2 et 713 respectivement, cig 7 ^ 0 est un nombre 
réel et £ est un petit paramètre. Nous montrons aussi dans ce chapitre que le résultat 
obtenu dans notre travail [41] est une meilleure estimation que le résultat de [30]. A 
la fin du chapitre nous donnons quelques exemples de comparaison entre les deux 
résultats. Nos résultats sont les suivants. 

Théorème 4.1.1 Pour |e| > 0 suffisamment petit, le nombre maximum de cycles 
limites du système différentiel (4-%) d u ' L bifurquent d’un centre linéaire x = — y , y = x 
en utilisant la théorie de moyennisation du premier ordre est 

Ai = max 1[ : y! 1 j • 

Théorème 4.1.2 Pour |er| > 0 suffisamment petit, le nombre maximum de cycles 
limites du système différentiel (f.2) qui bifurquent d’un centre linéaire x = —y, y = x 
en idilisant la théorie de moyennisation du deuxième ordre est 

A 2 = max {[H] + [açl] , [m] + [ç] - 1, [ail] , pys] , [»] + [ç] 

im + m 4¥i. pç*]. [“Pi+[= ÿi]+ m, 1 } , 

où K = min { [ÜÇ1] ,[¥]}• 

4.2 Preuves des résultats 
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Prenons 6 comme une nouvelle variable indépendante, le système (4.4) devient 

/ m \ 
dr 


d8 


= e cos#X e * rî cos* 9 — sin9P(9, r) + 0(e 2 ) 


i =0 


= eF 1 (0,r) + O(£ 2 ). 
Calculons maintenant i ? io(r) 


2 tt 


1 


.Fio(r) = I ( cos$X e 4 r * cos* 9 — sin9P(9, r) J dO 
n V 2=0 / 


1 m 

i y e , 

2 v 


2?r 

„4 f „„„!+1 , 


r* / cos* +i dd#-— Xair* / cos* 0 sinddd- —S^bjr 

2ir^ / 2 tt^ 


441 /. 442 


i =0 


,4+1 


0 


4=0 


0 


4=0 


2tt n ^ 2 tt 2tt 

x [ cos* d sin 2 9d9 — f cos* 9 sin 3 9d9 - ^d^r 3 f sm 4 8d8. 


2tt‘ 


o 


4=0 


2 tt 


0 


0 


Maintenant, en utilisant les expressions des intégrales en annexe A (notez que 
«fc+i = (2k + 1 )«fc), nous obtenons 


[^] 


*io(r) = 771 


2 tt 


X] e 2 *+l r 


2 tt 

2i+1 /cos 2 * +2 


4=0 


[^] 2 -rr 

8d9 - -ÿfcir 2,+1 [ cos 2 * 9 sin 2 ddd 

27 r t^ J 


--d 3 r 3 


[^] 


= r 


E a 2+1 2 i 

—e 2i+1 r 


[¥] 


4=0 


• „2i ^ 1 3 


4=0 


2 *+! (* + !)! 


^ 2 P’ - ~dkr a . 


4 = 0 


2* +1 (i + 1 )! 8 


Donc 


[^1 




w = r Epiî>"*- r E 


a ’- ’ „2 4 3 1 3 


4=0 


4=0 


2 i+1 (i + 1 )! 8 


& 2 ir - ô d o r • ( 4 - 5 ) 


Alors le polynomiale i ? io(r) peut avoir au plus Ai = max {[^], 1} ra¬ 
cines positives. Par conséquent, le système différentiel (4.2) peut avoir au plus 
Ai = max{[^], [^], 1 } cycles limites. 


4.2.2 Preuve du théorème 4.1.2 et certains lemmes 

Pour démontrer le théorème 4.1.2, nous allons appliquer la méthode de moyenni- 
sation du seconde ordre au système (4.2). Si l’on écrit / 1 (x), g 1 (x), h 1 (x) et l 1 (x) 
comme en (4.3), et 

ni n 2 n 3 m 

f 2 (x) = 5>x*, 5 2 (æ) = h2 ( æ ) = X s * æî ’ ( æ ) = 

2=0 2=0 2=0 2=0 
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Le système (4.2) s’écrit en coordonnées polaires (r, 9), r > 0 sous la forme suivante : 


r = e(cos 9 e^r* cos* 6 — sin 9P(9, r)) + e 2 cos 9v t r' 1 cos* 9 — sin 9Q(9, r ) , 


9 = 1 — f: cos 9P(9, r ) + sin ^y^r* cos* 9 — (cos 0Q(0, r) + sin 0 


<2 ( 0 , r) = X>* cos* 0 ^ 


gir * +1 cos* 0 sin 0 + yy r * +2 cos* 0 sin 2 0 + d^r 3 sin 3 


J + 2 fl ni., 2 fl L „• 3 , 


Considérons 9 comme une nouvelle variable indépendante, alors le système (4.6) 
devient 


— = e (cos 9eiP cos* 9 — sin 9P(6, r)) + e 2 ( cos é?yy r* cos* 9 — sin 9Q(9, r ) 

V i =o 


eiP cos 1 9 — P 2 (9,r) +P(9,r) 


(2 cos 2 9 — l) 


cos 9 sin 9 


x^Je^r* cos* 9 J + 0(e 3 ) 

i=0 ) 

= eF 1 (9,r) + £ 2 F 2 {9,r) + 0(£ 3 ). (4.7) 

Afin d’appliquer la théorie de moyennisation du second ordre, Fiq doit être identi¬ 
quement nulle. On a Fio = 0 si et seulement si 

b 2i = (2 i + 1) e 2 i+i, if 0 < i < fi and i / 1, 

< b 2 i = e 2 i+i = 0, if fi + 1 < i < Ai and i / 1, 
b 2 = 3e3 — 3dQ, if i = 1. 


ou fi = mm 


im.if]} 


Nous avons 


dFi(6,r) 

dr 


^ "i e i rl 


y^ia^r* 1 cos* 9 sin 9 — y^ (i + 1 ) 6 jr* 


x cos* 9 sin 2 9 — y^ (i + 2) Cjr * +1 cos* 9 sin 3 9 — 3 d^r 2 sin 4 9 
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ni 


n 2 — 1 
2 


^ 2 ie 2 «r 2 * 1 cos 2 * +1 0 + 3 e 3 r 2 cos 4 0 — y+a+ 4 cos* 0 sin 0 — ^ (2i + 2 ) 

2=0 2=0 2=0 

n 3 

x 6 2 i+ir 2?+1 cos 2 * +1 9 sin 2 6 - 3 b 2 r 2 cos 2 9 sin 2 9 - ^ (i + 2) c ++ 1 cos* 0 sin 3 9 

i= 0 
M 

—3do ?’ 2 sin 4 0 — ^ 6 2 jr 2 * cos 2 * 0 (( 2 i + 1 ) sin 2 0 — cos 2 0 ) , 


i=0 


É» 


et yi{9,r) 


jF\ (t, r) dt. Pour ce faire, nous réécrivons 
o 


[ m — 11 1" m ' 

L 2 J L 2 . 

22+1 „ 2 i +2 , 


Fl ((9, r) = ^ e 2 i+ir 2 * +1 cos 2i+2 9 + ^e 2i r 2 * cos 2i+1 0 + e 3 r 3 cos 4 ( 


i=0 


i=0 


n- 2- 1 

2 


-^a*r* cos* 9 sin 0 — ^ 6 2 j + ir 2 * +2 cos 2 *+ 4 0 sin 2 9 — dgr 3 sin 4 9 


i =o 


i =0 


L 2 J n 3 

-^ 6 2 jr 2 * +1 cos 2 * 9 sin 2 9 — & 2 r 3 cos 2 9 sin 2 9 — ^c ++ 2 cos* 9 sin 3 9 
2=0 2=0 


ni 


++ 2 t r 2 * cos 2 * +1 0 + e 3 r 3 cos 4 9 — ^a+ cos* 0 sin 0 
2=0 2=0 


2 


n 3 


X] 6 2 i+ir 2 * 


2î+2 „„ 0 2*+l 


0 sin 2 0 -^ 


i =0 


' cos"' 1 * V sin" y — ^c + +2 cos* 0 sin 3 0 — dgi ’ 3 sin 11 
i =o 

M n- 


M 


- 6 2 r 3 cos 2 0 sin 2 0 - ^ 6 2 ;r 2 * +1 ™° 2 * - 1 - ^ + 2 .,„ 2 *+i ^ 0 2 i +2 


^E4^T^+-cos-0. 
î= 0 ^ 2z + 1 


cos 


En tenant compte du fait que 

1/2 i\ 2i + 2 1 / 2i + 2 

22ï V * / ~~ 2ï + 1 2 2i + 2 V i + 1 

en utilisant les intégrales de l’annexe et comme 6 2 = 3e 3 


0 = 0 , 

- 3dg nous obtenons 


yi( 0 ,r) = ^+ 2 i r 2 * ^ 7 ^ sin( 2 Z + 1)0 +-e 3 r 3 ( 2 sin 20 + sin40) 

i =0 1=0 


n 2~l 

2 


'‘■l L “ J -i 

■^+ r ‘ (1 - cos<+1<, ) + E irr^i+ir 

i =0 i =0 


21+2 cos 2 *+ 2 0 sin 0 
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n 2 ~ 1 
2 


^ (0i\^ b2i+ir2i+2 S m sin ( 21 + V 0 ~ 5> 


n 3 


,i+2 


î=0 


(2z + 3) 


1 


1=0 


i =0 


1 


x cos î+1 0 + -—- cos* +3 I 
i - 1 - 3 

x sin( 2 Z 0 ), 


M i+l 

+ ^4r 3 (2 sin 29 - sin 40) - b 2i r 2i+1 ^ (3 iy 


(i + l)(i + 3) i + l 

i+l 

D 

l 


i =0 


/=1 


OU 


À,/ = ' 

Pour déterminer la fonction 


Pi,l - 


( 2 i+l) 
2 (i+l)/?j+i,i+i 
( 2 i+l) 


-, 0 <l<i, 

l = i + 1 . 


2n 

F 2 o( r ) = ^J dFl ^ ,r) yi & r ^ d9 + ^J W, r W 


2 tt 


notant par 


J ( r ) = 


0 


27T 

1 f fdFifrr) 


2ir J \ dr 
0 


2 ?r 

2/1 (<9, dé», I/(r) = ^ J F 2 (0, r)dd. 


Maintenant, divisons le calcul de la fonction I(r) en six morceaux. Pour cela, nous 
posons 

I(r) = h(r) + I 2 {r) + h{r) + h(r) + h{r) + I G (r), 


ou 


2 w m 

-1 Æ. 


h(r) = I ^ ? ’ a i r * 1 cos* 9sin9yi (9,r) d9, 
71 i i =o 


'*> = t 


2 -k 


/hF 1 


f ^ (2i + 2) & 2 i+i »" 2 * +1 cos 2î+1 0 sin 2 9 + 3b 2 r 2 cos 2 0 sin 2 y\ (9, r) d9 , 

n i =0 


-1 

/ 3 (r) = — (i + 2) ar l+1 cos 1 9 sin 3 9-yi (9, r) d9, 

71 o i=0 

2tt 

I 4 (r) = J3dlr 2 sin 4 9yi (9, r)d9, 


o 

2 ? / [f ] 


4M = s 


S ^2ie 2 ir 21 1 cos 2î+1 9 + 3 e 3 r 2 cos 4 9 | yi (0, r) d9, 
i =o 

/fj(r) = ~ô~ f ^y^& 2 ii’ 2a cos 2i 6 * (( 2 i + 1 ) sin 2 9 — cos 2 9)^j y\ (9 , r) d9. 
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Les intégrales seront calculées dans les lemmes suivants : 


Lemme 4.2.1 L’intégrale I\(r) est un polynôme de la variable r donnée par 

» 


— 1 
2 


J l( r ) - ^ A i,j a 2 ib 2 j+ir 2l+2j + 1 ~ ^2 niL^ + , L, fl 2 i +1 


i=0 j =0 


1=0 

"lirml 
2 J [ 2 J J 


2 *+ 3 (i + 3) ! 


x(3ai+ido + ai+2e3)r 2î+3 — EDE lj,iCi,ia2je2ir 2l+2 ^ 1 

i=0j=0 1=0 


1 ' 2 J /* / 9 - , __ 

+ J2 O zZ4;ê , M Q2»+l&2 J T 2i+2j+1 , 


7=0 7=0 


J =1 


oti Ajj est constante non nulle. 


Preuve. En utilisant les intégrales de l’annexe, nous obtenons 
{ a k +2 = (2/c + 3) (2k + l)ufc) 

2n 

( U1 \ I U1 ^ 

4+ 1 0)1=0. 


î 

(ai) / ( cos* 0 sin 0 J 1 ^ ^ r J ( 1 

\l=o J + 


4=0 


2 tt 


— COS 


(bi) f ( 1 cos* 0 sin( 




4=0 


V 


V —— 6 2 j+ir 2j+2 (- cos 2 ^+ 2 6 sin 0 
l=o 2j + 3 


;[ r V 


+y^ 7 j,i sin( 2 / + 1)9 j J = ^ ^ A^a2*fr2j+rr 2 * +2j+l , 

Z=0 / ) i=0 j=0 


OÙ .Aj j = 


2 tt 


2 j + 3 


(cr) /'(!>•' 1 

n Vi=0 


É7«Cy 


a î+l + l 


4=0 


2 *+l+ 1 (i + j+ 2 )! 


cos u sm ( 


' n 3 

E 

4=0 


c,r 


i+2 


(î + 1) (i + 3) 


cos ?+1 6 + . ^ „ cos î+3 0 ) ) = 0 . 


i + 1 


i + 3 




o 

ni —1 


/ 1 cos* 0 sin 0 ^ ( 2 sin 20 — sin4é*)^ 


3(2z + l)«i+l ,1 2i+3 

Z__, O i+0tA I QU a 2ï+l«0 r 


1=0 


2*+ 3 (i + 3)! 


(4.8) 
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2 tt 


( e l) 


2 vr 


ni \ / [ 2 ] j 

C ° S * ^ Sal ^ I X^^X 7 ^ 11 ^ + 1 )° 


i =0 


V 


3=0 1=0 


mm j 


- y^y^y^J j,lCi,ia2ie2ir 2l+2j i . 

i =o j =o 1=0 

27T / S 


(fi) J ^r\air l 1 cos* 0 sin 0 j (|e 3?’ 3 ( 2 sin 20 + sin40)) 
(2 i + l)ai +2 


ni — 1 
2 

E 

3=0 

27T / 


2*+ 3 (f+ 3)! 


®2i+1 e,3'r 


2i+3 


ni 


3 +1 


(gl) / ( ^X* r * 1 cosi $ s i n ^ X & 2 * r 2 j+ 1 X^ 7sin ( 2Zér) 


o 


vî=0 


0 =° 


i=i 


ni — 1 

““2“ fl 


■ * J ^ /rj • i i \ j + 1 

XI X—2—ÈÂ-Ai‘** + 1 ^^ +1 . 

i =0 i=0 Z=1 

On sait que la somme des intégrales (ai), ( 51 ) est le polynôme (4.8). 

Lemme 4.2.2 L’intégrale I 2 (r) est un polynôme de la variable r donnée par 

M"# 1 ] 


/ 2 (r) = X X 

i =0 j =0 

Tni — il 


(j + l)a i+j+l 


£ 


2*+i+i (2z + l)(z + j + 2)! 

3a i+2 


a 2 ib 2 j + ir 2 ‘ l+2j+1 


2 i + 4 (f + l)(i + 3)! 

n 3 


2=0 
ni 


+g 8(ïTï)“ iter ‘ +2 + Eït; 


a2i+i^2î’ 2î+3 

3 


^ 4 (z + l)(z + 3 ) 




i+4 


i=0 


îdff] 

S §2 i +i+ 2 (2j + 3)(2j + l)(i+j + 3)! 


(f + l)(10j + 4z + 15)a i+ j +1 


62 ,:+ic 2 i r 2i+2j+3 


V _ 3(5i + ll )«. t+2 _ 2i+5 

^ ?i +5 U 4-1 i (i + 9 \(i + 4U W2L2i + 1 


f^ 0 2 i+5 (i + l)(z + 2)(z + 4)! 

Preuve. En utilisant les intégrales de l’annexe, nous obtenons 
(a fc+ 3 = (2 k + 5)(2 k + 3) (2k + 1 )a k ) 

2 TT 

1 


( a 2) ^ 


^ J |x 0 + 1 ) cosJ $ sin 2 (l — cos 2i+1 0 )^ 


(4.9) 
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«1 q L 2 J q 

= g^ITT)^’'‘ +2 - g ¥^fwTsÿ awb2r2W 
mM ,.^ n 

-g g 


J= 

2 îr 


(b 2 ) / J Z (j + 1 ) bjri cos J 0 sin 2 0 
o V=° 


2ÏT3 6 2mr 2i+2 : 

V 


i =0 


cos 2 * +2 0sin0 + sin(2Z + l)0^j ^ = 0. 

i 7( n2 \ ( m . / 

(c 2 ) — ( z (* +1) œs * 0 sin2 0 ) Z C2 ^ 2j+2 ( p 

n Vr 0 / 7=0 


V 


2 j + 1 


cos 2j+1 0 H--— cos 2j+3 0 ) ) = Z TT 

2.7 + 3 // 4(i 


i=o 

«3 


^4(i + l)(i + 3) 

^ijrnâi 

\ ^ + 1 )( 1 °j + 4z + jjgb+2+l h r 2i+2j+3 

Z. Z. 2 <+j +2 {2 j + 3) (2J + !)(?■ + j + 3)! ° 2i+lC2j 


(j + 1)(j + 3) 

& 2 Cj ?’* +4 


n 3 ~ 1 


3(5? + ll)qj + 2 , 2 i +5 

Z 2 i + 5 (i + l)(i + 2)(i + 4 )! &2 2î+i 


7=0 


(d 2 ) — [ { — Z (* + 1 ) hr 1 cos* 0 sin 2 0 ) Zdor 3 (2 sin 20 — sin40)^) = 0 . 

21 i \ U / Vs > 


2 n 


M / -£ (i + 1 ) bir z cos* 0 sin 2 


7=0 


([f] J 

Z^’Z™ sin ( 2 ^ + 

^=0 z=o 


= 0 . 


(f 2 ) Z gg (j + 1 ) bjr * cos* 0 sin 2 0 ^j ^e 3 r 3 (2 sin 20 + sin 40)^ = 0 . 


r n2 \ / M \ 

(ë 2 ) / Z^' + 1 ) h i r3 cos J 0 sin 2 0 I yg 6 2 7 r 2 * +i ygÂ, ; sin( 2 Z 0 ) J = 0 . 

0 \J=0 / \i=0 1=1 J 

On sait que la somme des intégrales (a 2 ), ( 52 ) est le polynôme (4.9). 
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Lemme 4.2.3 L’intégrale I${r ) est un polynôme de la variable r donnée par 


mi?] 


h(r) = -EE(* + 1 )E lj,lKi,i c 2ie 2 jr 2l+2j+1 

i =0 j =0 1=0 

r w 2"i j p i rntiii 

+ E 22Bijb 2 i + iC2jr 2z+2j+3 - E ^i+5(f^4)T( 7d o + ( 4 * + 9 ) e3 ) 


i =0 j =0 


î=0 


w i+5 +E E 


n.3 — 1 

M L ~ J (2j + 3)* +1 




,2i+2j+3 


(4.10) 


i=0 j= 0 


{=1 


oti Bi j est constante non nulle. 


Preuve. En utilisant les intégrales de l’annexe, nous obtenons 

2tt 


«3 


ni 


( & 3) / Eü' + 2 ) 


ri+3 rnsi 


c ? r-' 1 cos J y sm 


in 3 9 j —— f l (l — cos î+1 0) = 0. 


<3=0 


1 //*• 

( fa 3) Y J Eü+») 

0 V=° 


.ts i+1 


Cjr JI+1 cos J 0 sin 3 0 


\ (l^l , 
E 2 Æ 3 

/ V *=° 


&2i+ir- 2i+2 x 


— cos 2î+2 0 sin 9 + 22^1,1 sin(2 1 + 1)9 




1=0 


Y Y Bi ’j b 2 i+l C 2 jr 2t+2j+3 - 

i=0 j= 0 


U + 1) 


où Bj i = 

1,3 2i + 3 


Y^ iK „ 


3a 


3,1 


i+j+l 


4=0 


2*+i+ 2 (i + j+ 3)! / ' 


2tt 


«3 


«•3 


( c a) / Eo'+ 2 ) 


ri +1 rnsi 


a=o 


Cir J 1 cos J sm 


1 


in 3 0 ( E c * r 


î+2 


4=0 


(i + l)(î + 3) 


-1 -cos ,+1 9 + 1 -cos î+3 9 ) 1=0. 

i 1 ?' “t - 3 


2tt 


(ds) j ^ ^ (z + 2) ar l+1 cos* 9 sin 3 9 J ^^dg?’ 3 (2 sin 20 — sin 40)^ 


0 

2 

E 

i=o 

2?r 


21aj+2 
2 i+5 (i + 4)! 


C2i+ldo r 


l„2i+5 


1 7 / ™3 \ f[f\ 3 

(es) 2^/ "E (i + 2) Cir l+1 cos 1 9 sin 3 6 ) E e2 i r ^E 7 il sm ( 21 + 1 ) 6> 
n V i =0 / i=o i=0 
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££< i + 1) £ 7 ^, 0 *^ 


, r 2i+2j+l 


i =0 j =0 


(fs) J (i + 2 ) Cjr l+1 cos* 6 sin 3 o\ ^e 3 r 3 ( 2 sin 20 + sin40)^ 


3(4 i + 9)«j + 2 2 i +3 

2 ,:+5 (i + 4)! 2 * +1 3 ' 


i 2 U n 3 

fe) ^ Eü' + 2 ) 

0 V=° 


Cir - 7 " 1 " 1 cos J 0 sin 3 i 


0 XV 2 ?+ 1 E/Mn( 2 ^) 


M L 2 J /rj • i q\ ®+l 

£ £ ^£ 

i =0 j =0 /=! 


^Â,/^j ) / 62 iC 2 jr+ir 2î+2j+3 . 


On sait que la somme des intégrales ( 03 ), ( 53 ) est le polynôme (4.10). 
Lemme 4.2.4 L’intégrale Ia{v) est un polynôme de la variable r donnée par 


I±{r) = — Y^ -r ———^7 -— a 2 j+i(inr 2î+3 

v y f^ Q 2 î+4 (i + l)(i + 3)! + 0 

77.1 q 77.3 q 

+ Ë»ô , p ai ^ !+2 + ëiTTTWTl 


^ 8 ^ + !) 


^4(i + l)(i + 3) 


Cjdor * +4 


9(7 i + 13)a i+ i 
2*+ 5 (z + l)(i + 2)(i + 4)! 


c 2 i+ idor 2î+5 , 


(4.11) 


Preuve. En utilisant les intégrales de l’annexe, nous obtenons 


( & 4 ) ^J (3 dy sin 4 6 ) i (! “ cos * 6 ) = Y. 


; 8 (* + i) 


aidr,r l+2 


9«i+l 2 

2 i + 4 (i + l)(i + 3 )!° 2i+1 ° r 


( 64 ) f (3dgr 2 sin 4 9) Y"' —- b 2 i+ir 2l+2 [— cos 2t+2 8 sin 6 

/ ZiX \ O 

^ „•_ r\ 


+y^jjj sin( 21 + 1)0 = 0 . 
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27T / n ^ 

( C 4) 2 n / ( 3ci o r2 sin4 e ) E^ î+2 


+=0 


—-—-- — . — cos* +1 0 + . — cos ?+3 0 

(i + l)(i + 3) i + 1 i + 3 


«3 


rcS" 1 

2 




^4(i + l)(i + 3) 


Cj4r î+4 - E 


9(7i + 13)aj+i 


i =o 


2 i + 5 (i + l)(i + 2)(i + 4) 


-C2î+l(fo 7 


,2i+5 


(d 4 ) 


2tt 

^ y (—3dgr 2 sin 4 0) ^Eg?’ 3 (2 sin 20 — sin 40)^ = 0. 


o 

2?r 


/[f 


( e 4 ) ^ f (~3dlr 2 sin 4 0) E^^Èm sin ( 2 ^ + 1 ) é * = °' 

'J • _ r\ » c\ 


o 

2tt 


i=0 Z=0 


(fi) E y (—3dgr 2 sin 4 0) ^e 3 r 3 (2 sin 20 + sin40)^j = 0. 


2îr 


i+1 


(g 4 ) / (3dg?’ 2 sin 4 0) I E &2 * r2î+1 E^4 sin ( 2 ^) ) = 0. 

o V *=o z=i / 

On sait que la somme des intégrales (a 4 ), ( g 4 ) est le polynôme (4.11). 
Lemme 4.2.5 L’intégrale I§{r) est un polynôme de la variable r donnée par 


Mf] 


= EE 


JQ+t-j+i 


î=o j^2^(2i + l)(z + j + l)! 


d2ie2jr 


2i+2j-l 


n 1 

-E 


9 


^8 (i + 1) 


mm 


aje 3 r 


1+2 


+ E 


3< ++3 2i+3 

-«27+1+3+ 


^ 2*+ 4 (i + l)(i + 3)! 


j(6i + 4j + 9)a i+j +i 

J=0 ^ 2i+i+1 (2i + 3) (2i + 1 ) (i + j + 2) ! 


+££; 


C2i&2j 


„2i+2j+l 


+ £ 


9( ‘ + 3)a,+3 - C2i+iesr zi+5 


2* +5 (i + l)(i + 2)(i + 4)! 


«•3 


9 


i +4 


'^4(i + l)(i + 3) Cie3r ' 


Preuve. En utilisant les intégrales de l’annexe, nous obtenons 

(ofc+ 4 = (2i + 7) (2i + 5)(2fc + 3)(2fc + l)a fc ). 

2t r 

<«> è 


(4.12) 


ni 

E j e j r ^~ 1 cos J+1 0 J | — E E -, + a (40 
O=o / V i=o l + 


— cos 


2i+1 0 ) 
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n i 

-£ 

i =0 


9 


8 (* + l) 


Oie 3 r î+2 + ^ 

i =o 


3u»+3 2/ ( 3 

-j-<22i+ie3r 2 +,i 


2 i+4 (i + l)(i + 3)! 


¥][?] 


+ EE: 


J " î+J+1 a 2 i e 2 j r 2 *+ 2 ^ 1 . 


i= 0 ^ 2 i+j ( 2 * + 1 )(® + i + 1 ) ! 


1 11- \([^] x 

( b s) ^ / E-^ rJ 1 cos:,+10 - E wr 2î+2 [ 

0 V=° / V i=0 


— cos 2î+2 0 sin 0 


+^ 7 i,l sin( 21 + 1 ) 0 ] = 0 . 


1=0 


<C5> 2 Ï 


k 7 

0 V'=° 


cos- 


( 1 + 1 , 


-E c * r 


n 3 


i=0 


i+2 V - 

{{^(i+m+sr 


c;r l+2 


i + 1 


mm 


cos * +1 0 H-cos î+3 0 ) ) = — — 

i + 3 II ^Ad 


«-3 


9 


i =0 


4(i + l)(z + 3) 


Cie 3 r 


i+4 


lW j( 6 z + 4j + 9)a j+ j + i 21 + 2 . 7-1 

^ ) f^2 i + j + 1 (2i + 3)(2i + l)(i + j + 2)\ 21 2] 


2 


\ ^ 9 (i + 3)a i+3 2i +5 

+ E 2 i +5(i + l)(i + 2 )(i + 4)! C2i+ie3? ' 


<*>2Ï 


1 cos î+1 0 ^j ^^o r3 (2 sin 20 — sin 40)^ = 0 . 


o 

2 ?r 


([*: 


(e 5 ) ^ / f E je l ri 1 cosi+1 0 É e 2 ir 2 *Z%l sin( 2 1 + 1)0 = 0 . 


o 

2 ?r 


0 =° 


i=0 1=0 


(fs) ^y" 1 cos î+ 1 0 ^ ^e 3 r 3 (2 sin 20 + sin 40)^ = 0 . 


2 tt 


i+1 


te) / E je i rJ lcosJ+l0 -E^^Eâ^w) = °- 


0 =° 


i=0 Z=1 


On sait que la somme des intégrales (as), (g 5 ) est le polynôme (4.12). 




64 


Chapitre 4. Nombre maximal de cycles limites pour les systèmes 
différentiels polynomiaux généralisés de Kukles 


Lemme 4.2.6 L’intégrale I&(r) est un polynôme de la variable r donnée 

7.=0 1=0 v J ' 


par 


i=0 j=0 

r ”3-! 
U l 2 


+E E 

i=0 j=0 

Preuve. En utilisant les intégrales de 


3Ui+j+l u oi+2j+3 

2 i+i+3(i + j + 3)\ b2l2j+1 


(4.13) 


annexe, nous obtenons 


(<*) è f E b 2 jr 2 j cos 2 - 7 0 ((2j + 1) sin 2 6 — cos 2 0) J | . 0,1 r 1 (40 

o \i=o / Vi=o* + 


— cos* +1 


2W ^ + 2)! ^+^+ 


■»))= E E 

î=0 j=0 

2n ( U 

(bg) f f ^ 02 jî’ 2 - 7 cos 2 - 7 0 ((2j + 1) sin 2 0 — cos 2 0) 

w7r, Z \ 1=0 


x 


o V=° 

Z 

— cos 2î+2 0 sin 0 + ^7i,z sin(2/ + 1)0 

2ir 


/[^ , 
E 


v 


i=0 


2i + 3 


02i:+ir 2i+2 


1=0 


= 0. 


(cg) fj jr 21 cos 2t 0 ((2i + 1) sin 2 0 — cos 2 0)^ ^ 

x ( ——— — c 1 r 7+2 —t — r cos j+1 6 + -— - cos- 7 " 1 " 3 0^ ^ 


4+2 


t--t -rCo-r-' 1 -cos J a -\ - 

0 + 1)0 + 3 ) J + 1 J + 3 


u [ 2 J o 

V V _ ctot+ - 7+1 _5 9îC9 , , , r 2 *+ E + 3 

Z. Z. 2*+i+ 3 (f + j + 3)!^ 2 C2j+1 

2tt 


(ds) M {-P** cos 2 * 0 ((2z + 1) sin 2 0 — cos 2 0)^ ^dg?’ 3 (8 sin 20 


— sin 40)) = 0 


2tt / ^ \ 

(eg) j ( — ^ 02 i^ 2 * cos 2 * 0 ((2i + 1) sin 2 0 — cos 2 0) 


o v 7=0 


/[f] 

Ÿ,e 2j r 2j 
\j =o 


sin(2/ + 1)0^ = 0. 


1=0 


(fg) J S ^p 2 jr 21 cos 2 * 0 ((2i + 1) sin 2 0 — cos 2 0)^ ^e 3 r 3 ( 8 sin 20 


o v 7=0 
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+ sin4#)) = 0. 

(gfj) f ('y^b 2 ir 21 cos 2t 6 (( 2 i + 1 ) sin 2 6 — cos 2 6 ) j | ^ 62 jr 2jl+1 

0 V'=o J \j =0 

i+i 

sin(2Z0) 
z=i 

On sait que la somme des intégrales (ag), (gg) est le polynôme (4.13). 
En utilisant les lemmes 4.2.1-4.2.6, on obtient : 



2tt 

/(r) = èri ( ' dFl dr r) yi (0 ’ r) ) d6 = /l(r) + h{r) + /s(r) + /4(r) + h{r) + /6(r) ‘ 

0 


Maintenant, nous considérons l’intégrale II{r ) 
3dQ, d’ où II(r ) est égal à 


2tt 


JF -2 (?’, 6») dd. On a & 2 


3e 3 - 


E 


CKi+l 


î“2i+i r2i+1 -Ea 


®.r 2i+1 - Idgr 3 


2 ‘ +1 (ï + 1)! =-2‘+ 1 (i+l)! 

, Mm, , 

_y^ «i+i+lQ 2 i 02 j+l r 2 i+ 2 . 7 +l \ ^ \- -io, . ; , ; tf‘2i ■ l'g, T , 2 »+ 2 j +3 


i=0 j=o 
—1 


2 i +i+ 1 (i + j + 2)\ 


E E 2 i+j+2(i + i + 3)!' 


+ E 


Oti +1 


i=0 

"3-1 


2* +2 (i + 3)! 


a 2 i+i ((* + 1)4 + (* + 5)(2i + 3)e 3 ) r 


,2i+3 


mm 


+ E 


ati+2 (4 - (* - s ) e 3) C 2 i+l ^ 2i+5 1 y-^ (j + ./)<>, . //2/f 2,/ 2 i+ 2 j — 1 


i=0 


fflî 


+EE 

i=oj=o 

n-^ —1 


2 i+3 (i + 4)! 


(?■ + j l)aj + jC2ie2j 2i+2j+l 


+ EEï 


2 *+l+ 1 (i + j+ 2 )! 


i=01=0 

r n 3 — 1 

m L 2 

EE 

l=o 1=0 


2*+l(i + j + 1)! 


(5i j + 3)aj+j+ib2iC2j+l r 2i+2j+3 
2 *+i + 2 (2i + 1 ) (ï + j + 3) ! 


2 j M 

+ E E 


(* j) a i+j+1 


^ ^2^+ 1 ( 2 j + 1 )(z+j + 2 )! 


a2i+i&2ir 2i+2 l +3 . 


(4.14) 
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De finalement, F 2 q (0 = I(r) + II(r) est égal à 


mm 


E E E hj a 2 i b 2 j+l r 2i+ ^ +1 + EE EijCl2i€'2j'ï' ^ 

z=0 j =0 j= 0 j =0 


— 1 
2 


3ai+ia2i+ido 2i+3 a ?+i (3do + 2(22 + 3)es) C21+1 

2 i+2 (i + 2 )! r ^ 2 *+ 3 (i + 3 )! 

Mi?] [?][?] 

+ E E ^,^ 2 *+ic 2 j r 2î + 2 ^ 3 + ^2^2F id C2ie 2 jr 2l+2j+1 

i =0 j =0 i =0 j =0 

+ E 


,2i+5 


«i+1 


[¥] 


î=0 
— 1 


2 i + 1 (?; + 1 )! 


^+ir 2i+1 - 


* Q2if 2i+1 ~ Ur 3 


i =0 


2 i+1 {i + 1)! 


1 

M L 2 


+ ^ ^^02^+162,r 2i+2j+1 + ]T H id b 2i C 2 j+ ir 2i+2 i +3 , ( 4 . 15 ) 
i =0 j =0 i =0 2=0 


, _ a _ (2z + j + 2) Qj+j+i 

' iJ “ 2 M+i ( 2 ?; + 1 ) (* + j + 2)! ’ 

, _ (2i 2 + 2j 2 + Aij + i + 2j)aj+j . 

“ 2 i +J(2* + l)(i + j + l)! * ho lj ' 1 ^ 

„ _ (4z 2 + 12j 2 + lOij + 19i + 34j + 24)aj_|_j_|_i 

i,j ~ i,j 2*+i+2( 2 j + 3) (2 j + l)(z + j + 3)! ’ 

^ (4z 3 + 8 j 3 + 20j 2 z + 22 j 2 + 16i 2 j + 32ij + 4f 2 + 12j — 5i — 3)oti + j 

i,j ~ 2 i +i+ 1 (2i + 3) (22 + l)(i + j + 2)! 

i 

- (i + DE^M. 

1=0 

r _ ( 2 ^ + 1 ) g ^ ,_ (2i + l)qj + j + i 

iJ ^ 2 Pi ’ Z M + 2 < +J+2( 2 j + 1)(* + j + 2)! ’ 

r _ (2j + 3) ~ ~ _ (4i - 2j + 3)aj + j + i 

i,j ^ 2 ^ j ' 1 2 < +i+3(2*+ !)(*+j+ 3)!’ 


Par conséquent, F20 a au plus max max { [^-] + , [^] + [y] — 

m + [t ] * m + m , [?] . [*p] + f 1} 


tives. 



racines posi- 
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4.3 Applications 

Dans cette section, nous montrons que dans les exemples suivants, la borne Ai 
(A 2 respectivement) est atteinte, (les calculs ont été vérifiés à l’aide de Maple). 

Exemple 4.3.1 (m = l,n 2 = 2 ,713 = 1 et m = 6 ) 

On considère le système 

x = —y + £l 1 (x)+£ 2 l 2 (x), 
ÿ = x - e(/ 1 (x) + g 1 (x)y + h 1 (x)y 2 + dly 3 ) 

-£ 2 {f 2 (x) + g 2 (x)y + h 2 [x)y 2 + d\y 3 ), (4.16) 

OÙ 

f 1 (x) = a 0 + aix, g 1 (x) = b 0 + b\x + b 2 x 2 , h 1 (x) = c 0 + c\x, 
l 1 (. x ) = eo + e\x + e 2 x 2 + e 3 x 3 + e 4 a ; 4 + e^x 5 + eçx 6 , 
f' 2 ( x ) = p 0 + pix, g 2 ( x ) = q 0 + qix + q 2 x 2 , h 2 ( x ) = s 0 + six 
and l 2 ( x ) = vq + v\x + v 2 x 2 + v 3 x 3 + u 4 a ; 4 + usa : 5 + vqx 6 . 


Tout d’abord, nous allons étudier les cycles limites de l’équation différentielle (4-16) 
en utilisant la théorie de moyennisation du premier ordre. De (4-5), nous obtenons 

F w (r) = ^(ei- b o+ 1 j ( b 2 ~ 3do + 3e 3 ) + ^e 5 r 4 ^ . (4.17) 


Maintenant, si nous fixons bo = 1, b 2 = 10, dç, = —4, ei = 3, e 3 = —4 et e 5 = -. 

Nous avons l ? io(r) = - (l — r 2 ) (4 — r 2 ) , qui a exactement deux racines positives. 

Dans [30], le système (4-16) lorsque l 1 = l 2 = 0 a exactement une seule racine 
positive. Afin d’appliquer la théorie de moyennisation du second ordre, nous avons 
besoin que Fiq soit identiquement nulle. Pins, à partir (4-17) nous considérons ei = 
b 0 i b 2 = 3e 3 — 3 d\, et es = 0. De Eq. (4-15), nous obtenons 

„ . , 1 , 1 1 2 (3 1 , 1 

F 2 o{r) = --Mo + a 0 e 2 - c 0 e 0 + -v\ - -qo + r I -d 0 «i ~ g&ic 0 - -c^ 


+ -e 4 a 0 + -v 3 - -d 0 - -q 2 +r —e 4 c 0 - -e 3 c 4 + — e 6 a 0 


l,i 5 A 25 6 

+ l6 d » C1 + Ï6 VS ) + 32 emr 


11 


15 


24 


Maintenant, si nous fixons ao = 0, ai 
dl = 2, e 0 = 0, e 4 = 3, e 3 = -4, e 4 


1, bi = 3, c 0 = -, ci = -1, dl = -4, 
5 8 

xx, e 6 = -xx, v\ = 0, u 3 = 1, v 5 = 2, 
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qo = —2 et q 2 = —. Nous avons que EWr) = — (l — r 2 ) (4 — r 2 ) (9 — r 2 ) , qui a 
45 36 

exactement trois racines positives. Dans [30], le système (4.16) lorsque l 1 = l 2 = 0 
a exactement deux racines positives. 


Exemple 4.3.2 (m = 2, n 2 = 3, n 3 = 4 and m = 7). 

On considère le système (4.2), oti 

/ 1 (x) = a 0 + aix + a 2 x 2 , g 1 (x) = &o + M + M 2 + M 3 , 
/i 1 (x) = co + ci» + c 2 x 2 + c 3 x 3 + C 4 X 4 , 

Z 1 (x) = eo + eix + e 2 x 2 + e 3 x 3 + e 4 X 4 + esx 5 + e 6 X 6 , 

/ 2 (x) = po + pix + p 2 x 2 , g 2 (x) = go + gix + q 2 x 2 + q 3 x 3 , 

h 2 (x) = so + six + s 2 x 2 + s 3 x 3 + S 4 X 4 

and l 2 (x) = xq + xix + x 2 x 2 + x 3 x 3 + X 4 X 4 + X 5 X 5 + vqx 6 . 


En utilisant la théorie de moyennisation du premier oindre. De l’équation (4.5) nous 
obtenons 

r ( 5 35 \ 

Fio(r) = ~ fei - 6 0 + r 2 (~b 2 - 3d[ + 3e 3 ) + -e 5 r 4 + —e 7 r 6 j . (4.18) 

116 

Maintenant, si nous fixons e\ = 3, 60 = 1, &2 = 10, dj = —4, e 3 = ——, 

Zj ( 

56 32 1 

es = — and e 7 =--. Nous avons F\o(r) = — (1 — ?’ 2 ) (4 — r 2 ) (9 — r 2 ) , qui 

45 315 36 

a exactement trois racines positives. Dans [30], ce système lorsque l 1 = l 2 = 0 a 

exactement une seule racine positive. 

Afin d’appliquer la théorie de moyennisation du second ordre, nous avons besoin 
que Fio soit identiquement nulle. Pins à partir de (4-5) nous considérons ei = 60 , 
& 2 = 3e 3 — 3 d,Q, es = 0 et e 7 = 0. De l’équation (4.15) nous obtenons 

_, . . 1, 1 1 2 (3 -, 1, 1 

F 2 o{r) = --Mo + « 0^2 - c 0 e 0 + -v\ - -qo + r I -d 0 ai ~ gMo - |Ciei 

3 1 3311 3 1 

P-e^ao - |C 2 e 0 + -v 3 - -d 3 - g <72 - -a 2 bi - -6300 + ^e 2 a 2 


1 


11 


1 


1 


7 


+ r To *" 2^1 + X7 e 4 c 0 - Txe 2 c 2 - xeic 3 - 7^&3°2 - Mo 


48 


24 


12 


48 


48 


5 1 1 15 1 5 \ 

+ Y^ e 4 a 2 ~ g C 4Co — g e 3 c l + “^“ e 6Û0 + Yg®0 c l + ) 


1 


35 


1 


+ r — e 4 c 2 - -—c 2 & 3 - —e 3 c 3 + —e 6 a 2 - —c 4 e 2 


64 


128 


64 


64 


25 3 f1 35 \ « / 7 

+ 32 ee<: ” + T28 C3 ‘ J » + i28” 7 J +r (sO^ 2 


16 

11 , 13 

-0 3 C4 — -C4e4 

1280 640 


+ 256 C4e6r 


10 
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25 135 1 

Maintenant, si nous fixons üq = 0, a\ = 1, ai = —, b\ = 3, 63 =-, co = -, c\ = 

54 118 5 

598719761 —472 

_1, C2 = -5, c 3 = toow , ccn > C 4 = —, dl = -4, d 2 = 2, eo = 0, e, = 3, e 2 = 1, 


63 = 


183274650 
116 1057617 

ee = 


2950 


27 ’ 226265 ’ 

—2 et qo = -. Nous avons que F 2 o(r) = 

675 


2673 

99 n 1 29818651106 

, vi = 0 , v 3 - 1 , v 7 - - 100 ^,-, 000 ^ r , qo - 


1 


14400 


12371038875’ 

(l — r 2 ) (4 — r 2 ) (9 — r 2 ) (16 — 


r 2 )(25 — r 2 ), qui a exactement cinq racines positives. Dans [30], ce système quand 
l 1 = l 2 = 0 a exactement deux racines positives. 




Conclusion et perspectives 


Cette thèse de doctorat porte sur un aspect important de la théorie qualitative 
des systèmes différentiels planaires, à savoir les cycles limites. L’importance de 
déterminer le nombre des cycles limites fait l’objet de la deuxième partie du 16 eme 
problème de Hilbert. 

La recherche des cycles limites des systèmes différentiels dépendant d’un petit 
paramètre peut être étudiée au moyen de la méthode de moyennisation. 

Nous continuons à travailler sur des problèmes analogues. On se propose 
d’étudier le nombre maximum de cycles limites du système 


I x = -y - l(x) - p(x)y , 

\ÿ = x - f(x) - g{x)y - h{x)y 2 - q{x)y 3 , 

où l(x) = el 1 (x) + £ 2 l' 2 (x), p(x) = ep 1 (x) + s 2 p 2 (x), f(x) = e/ 1 (x) + £ 2 f 2 (x), 
g(x) = £g 1 (x) + £ 2 g 2 (x), h(x ) = £h 1 (x) + £ 2 h 2 (x) et q(x) = £q 1 (x) + £ 2 q 2 (x). Pour 
chaque k = 1,2, l k (x), p k (x), f k (x), g k (x), h k (x) et q k (x), sont de degrés mj, 
m, ri2, ns et n 4 respectivement, e est un paramètre réel suppposé petit. 



Annexe 


Dans cette annexe, nous rappelons quelques formules qui seront utilisées au 
cours de cette thèse (voir pour plus de détails ([1], [20])). Pour i > 0 nous avons 

2t r 


J cos* 9 sin- 7 8d0 / 0, 
o 


si i et j pair, 

0, si i ou j impair, 


2tt 


cos* 9 sin- 7 9d9 = < 


7 

2 k ~ 1 kV 

TT a k 

2 k (k + 1)! 

37T a k 

2 k + 1 {k + 2)V 

\h-KOlk 


[ 2 k+2 {k + 3)!’ 
où ai = 1 • 3 • 5 • • • (2i — 1). 


si i = 2k et j = 0, 

, si i = 2k et j = 2, 

si i = 2k et j = 4, 

si i = 2k et j = 6, 


2tt 


2tt 


cos* é^sin^ 9 su\(2l9)d9 ^ 0, si i et j impair, 

T 

cos* #sin J #sin((2l + 1 )9)d9 / 0, si i pair et j impair, 

0, si i impair ou j pair, 


2tt 


cos* 9 sûP 9 sin((2Z + 1)9)d9 = < 


7T Cij, si i pair et j = 1, l > 0, 


nK h i, si i pair et j = 3, l > 0, 
où Cn e t Aj i sont des constantes réels non nul. 
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2 7T 


cos* 0 sin - 7 0 sin( 2 W)d9 = < 


0 , si i ou jpair 
7r(2 k + 1 )atk 


2 k {k + 2 )! 

irk(2 k + 1 )cnk 
2 k ~ 1 (k + 3)! 

3n(2k + 1 )ak 
2 k + 1 (k + 3)\ 

3ir{k — 1)(2 k + 1 )afc 


, si 7 = 2 fc + 1, j = 1 et l = 1, 

, si 7 = 2k + 1 , j = 1 et l = 2 , 

, si z = 2k + 1, j = 3 et l = 1, 

, si i = 2k + 1 j = 3 et l = 2, 


2 k (k + 4)! 

7i 'Cij, si i impair et j = 1, l > 0, 

^ 7 r Ki i, si i impair et j = 3, l > 0, 
où Cn and Ki,l sont des constantes réels non nul. 
e 

J cos* t s'm tdt = . ^ - (1 — cos * +1 9 ) , 
o 

e r. 

/ cos 2 * +1 tdt = ^ 7 ij sin( 2 Z + 1)0, (voir ([20], P153)). 

J i _ 


;=o 


cos 2 * +1 1 sin 2 tdt = 


cos 2i + 2 0 sin 0 + ^^ sin( 2 Z + 1 ) 0 , (voir ([ 20 ], 


P151)) 

6 


cos 2 * tdt = - 7 TT 


1 ( 2 i 


2 2 * \ i 


9 + sin(2Z0), (voir ( [20], P153)), 


i=i 


. . _ 1 / 2i + 1 \ 1 1 / 2i \ 1 

ou 2 2 * V i - l ) 21 + 1 ’ ^ 2 2 * V i + l ) V 

e 

/ 2 1 i 

cos* t sin 3 tdt = 7 - 7 - 7 -cos * +1 0 H-cos * +3 0, 

(i+ !)(* +3) i + 1 f + 3 


cos 2 1 sin 2 tdt = — (40 — sin 40) 

O Z 


sin 4 tdt = — (120 — 8 sin 20 + sin 40), 

O Z 
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e 

I cos 4 tdt = — (120 + 8 sin 20 + sin 40). 

J 32 

o 
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